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During  this  period  work  was  concentrated  in  the  fol  lowing  areas 


A.  “'Crystal  Growth 

B.  Hydrogen  /Kinetics  Theoretical^ 

... _ s^-  ..  —  -  ■ 

C.  ~r Hydrogen /Kinetics  r.* 

D.  '•.Devices^ 

F..  Review^  Paper 

F.  7>Far  IR  Spectroscopy^ 

G.  ’■Magnetic  Resonance  Study  of  defect  States, 

H.  -free  farrier  Optical  Absorption 

I.  Hydrogen  pi f fusion ? 

J.  'Electron  Trapping  States 


Experimental 


V 


) 


Each  area  will  be  discussed  separately  in  what  follows. 


A.  Crystal  Growth 


The  present  status  of  the  rutile  crystal  growth  project  is  extremely 
encouraging.  We  have  now  achieved  highly  reproducible  growth  rates  of  the 
order  of  1  mm/day.  The  resulting  crystals  appear  to  be  of  very  high  quality 
and  it  now  appears  almost  certain  that  this  growth  technique  will  indeed 
live  up  to  all  of  our  expectations  and  enable  us  to  produce  crystals  of 
mucli  higher  perfection  and  purity  than  have  ever  been  achieved  before.  A 
rather  surprising  number  of  problems,  some  of  a  very  fundamental  nature, 
had  to  be  overcome  to  r' ach  our  present  position.  Details  on  some  of  *hc 
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more  significant  problems  are  given  below: 

1 •  Transport  Mechanism 

Very  little  information  had  been  available  concerning  the  thermal 
stability  of  the  transport  agent  TeCl^,  which  is  obviously  of  great  importance 
in  understanding  the  transport  mechanism:  Firstly,  the  chemical  kinetics 
depend  on  the  availability  of  TeCl^  in  the  growth  ampoule.  Secondly,  a 
temperature -dependent,  reversible  dissociation  such  as 

TeCl4  $  TeCl2  ♦  Cl2  ,  (1) 

could  establish  density  gradients  in  addition  to  the  gradients  resulting 
from  the  temperature  profile,  and  thus  enhance  the  driving  force  for  con¬ 
vection  and  material  transfer.  It  was  felt  that  by  an  appropriate  choice 
of  the  growth  temperatures  (within  the  limits  irposed  by  chemical  and  inter¬ 
face  kinetics  considerations),  this  additional  transport  mechanism  could  be 
optimized,  if  a  more  quantitative  understanding  of  (1)  could  be  achieved. 

Thus,  calculations  of  the  kinetics  of  this  reaction  were  performed,  based 
on  available  thermo-chemical  data.  These  showed  that  thermal  dissociation 
according  to  (1)  indeed  occurs,  a  fact  only  suspected  hitherto.  The 
calcu’ition  was  checked  by  measurements  of  the  thermal  dissociation.  High 
purity  leCI^,  synthesized  in  situ  as  described  below,  was  heated  in  a  quart: 
bourdon  gauge  and  the  resulting  pressure  determined  over  the  temperature 
range  of  interest.  The  experimental  results  fit  the  theoretical  values 
very  well,  taking  into  account  the  limited  precision  of  the  available  heats 
of  formation  that  had  to  be  used  in  the  calculations. 

This  effort  has  given  us  major  insight  into  the  complex  transport 
mechanism.  Our  work  has  verified  that  the  transport  reaction  originally 
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assumed, 


TeCl4  +  Ti07  <-*•  Te02  +  TiCl4  , 

is  indeed  responsible  for  the  TiO,  transport.  However,  the  calculations 
above  show  that,  at  the  growth  temperatures,  the  TeCl^  is  more  than  90% 
dissociated  according  to  Eq.  (1).  It  is  interesting  that  the  loss  of  trans¬ 
port  agent,  due  to  thermal  dissociation,  is  more  than  offset  by  the  enhance¬ 
ment  of  convection  due  to  the  same  reaction.  This  clearly  demonstrates  the 
complexity  of  the  system  we  are  dealing  with  and  amply  justifies  our  con¬ 
tention  that  a  thorough,  detailed  understanding  is  absolutely  essential  to 
success  in  dealing  with  materials  like  Ti02> 

Various  transport  phenomena  which  we  had  observed  earlier  and  were 
unable  to  explain  with  mere  thermally  induced  density  gradients,  are 
qualitatively  understood  now.  Quantitative  studies  of  the  relative  weight 
of  thermally  and  concentration  induced  convection  are  planned.  Furthermore, 
with  little  more  effort,  we  will  be  able  to  contribute  thermodynamical 
information  about  the  TeCl^  system  with  substantially  higher  accuracy  than 
are  now  available.  Since  TeCl4  is  quite  frequently  employed  in  the  growth 
of  oxide  crystals,  this  information  is  of  general  importance. 

2 .  Purity  of  Transport  Agent 

A  major  experimental  difficulty  resulted  from  the  extreme  hygro- 
scopicity  of  TeCl4.  This  uncertainty  has  been  eliminated  by  synthesizing 
the  tellurium  tetrachloride  in  situ  from  high  purity  tellurium  and  high 
vacuum  distilled  chlorine.  This  refinement  has  dramatically  improved  the 
reproducibility  of  the  transport  rates. 

3 .  Profi le  Inside  Growt h  Ampoule 

Further  measurements  of  the  temperature  distribution  inside  the 
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growth  ampoules  as  well  as  in  simulated  open  geometries  led  us  to  a  major 
breakthrough  in  devising  an  appropriate  heating  geometry.  In  the  previously 
used  arrangement,  besides  growth  onto  the  seed,  heavy  deposition  of  TiO, 
occurred  also  on  the  wall  of  the  growth  ampoule  over  a  zone  of  several 
centimeters.  This  created  several  problems:  The  growth  rate  was  reduced 
by  reduction  of  the  degree  of  supersaturation  at  the  seed.  The  ultimate 
size  of  the  grown  crystal  was  limited.  Equally  importantly  for  present 
purposes,  the  deposits  on  the  ampoule  wall  obscured  visual  monitoring  of 
the  seed  growth.  A  novel  ampoule  geometry  has  been  developed,  in  which 
the  seed  is  cooled  by  radiation  to  somewhat  lower  temperatures  than  the 
quartz  wall.  The  technique  consists  essentially  of  weakly  coupling  the 
seed  crystal  to  a  low  temperature  bath  (room  temperature) ,  by  means  of 
light  pipe.  The  technique  relies  on  the  fact  that  the  thermal  emissivity 
of  the  TiO,,  is  much  higher  than  that  of  quartz.  Thus,  we  are  able  to 
produce  the  seemingly  paradoxical  situation  in  which  a  seed  crystal,  com¬ 
pletely  enclosed  in  the  quartz  ampoule,  is  maintained  at  a  temperature 
below  that  of  the  ampoule.  This,  of  course,  is  the  absolute  optimum  in 
growth  conditions.  Consequently,  growth  onto  the  ampoule  wall  is  com¬ 
pletely  excluded,  and  all  the  transported  material  is  deposited  on  the  seed. 
This  enables  us  to  substantially  reduce  the  seed-to-source  temperature 

l 

difference  and  hence,  the  temperature  oscillations  reported  earlier  can 
be  avoided,  even  in  a  vertical  arrangement. 

Two  papers  on  our  crystal  growth  efforts  were  accepted  for  presentation 
at  the  Second  National  Meeting  of  the  American  Association  for  Crystal  Growth 
in  August  at  Princeton.  The  abstracts  are  enclosed. 


-4- 


GROWTH  OF  Ti02  (RUTILE)  BY  CHEMICAL  VAPOR  TRANSPORT + 


F.  Rosenberger  and  M.  C.  DeLong 


Physics  Department,  University  of  Utah 
Salt  Lake  City,  Utah  84112 

Rutil.'.  single  crystals  of  typically  2  cm  diameter  and  2  mm  thick¬ 
ness  were  grown  homoepitaxially  in  closed  quartz  ampoules  employing  the 
overall  reaction.. 

Ti02  ♦  TeCl4  t  Te02  ♦  TiCl4  .  (1) 

Thermochemical  calculations  showing  that  TeCl,  dissociates  according 
to 

TeCl4  t  TeCl2  ♦  Clj  (2) 

were  confirmed  by  pressure-temperature  measurements.  The  role  of  (2) 
in  the  transport  mechanism  is  discussed.  Calculations  for  Cl2,  TeCl2 
and  TeCl4 :C12  mixtures  and  corresponding  transport  experiments  give 
strong  evidence  that  TeCl4  is  the  actual  transport  agent. 

The  growth  and  etching  form  of  rutile  obtained  in  (1)  was  found 
to  consist  of  {ill}  faces  --  in  disagreement  with  experimental  and  theo¬ 
retical  results  hitherto  reported. 

Preparational  procedures  are  described  and  details  on  the  growth 
process  which  was  performed  in  a  novel  heating-cooling  arrangement  are 
given. 


Work  sponsored  by  the  Advanced  Research  Projects  Agency  on  Grant 
No.  DAHC-15-71-G8. 

For  a  discussion  of  the  consequences  for  heat  transfer  in  CVT,  see 
paper  by  Rosenberger,  DeLong  and  Olson. 
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HEAT  TRANSFER  AND  TEMPERATURE  OSCILLATIONS 
IN  CHEMICAL  VAPOR  TRANSPORT^ 


F.  Rosenberger,  M.  C.  DeLong  and  J.  M.  Olson 

Physics  Department,  University  of  Utah 
Salt  Lake  City,  Utah  84112 

The  temperature  profile  in  a  crystal  growth  system  plays  a 
crucial  role  in  morphological  stability  considerations  and  interpre¬ 
tation  of  growth  phenomena  in  general.  In  growth  from  vapor,  particu¬ 
larly  in  closed  ampoules,  it  is  difficult  to  determine  the  actual 
temperature  conditions.  Thus,  the  temperature  distribution  of  the  heat 
source  is  commonly  assumed  to  also  apply  to  the  growth  system. 

Temperature  measurements  inside  a  closed  CVT  system  employing 
thermocouples  revealed  that  the  temperature  profiles  of  heat  source  and 
growth  ampoule  can  greatly  deviate.  Transport  agents  that  thermally 
dissociate,  such  as  TeCl^  and  can  lead  to  sizeable  heat  transfer 
rates  in  a  temperature  gradient.  Under  hydrodynamically  unstable  con¬ 
ditions,  temperature  oscillations  of  -  25^C  were  measured  in  the  gas 
phase. 

A  brief  discussion  of  heat  transfer  mechanisms  in  chemical  vapor 
transport  is  illustrated  with  measurements  under  different  geometrical 
and  compositional  conditions. 


Work  supported  by  the  Advanced  Research  Projects  Agency  on  Grant 
No.  DAHC- 15-71-G8 . 


See  paper  on  Ti02  transport  with  TeCl4  by  Rosenberger  and  DeLong. 
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B.  Hydrogen  Kinetics  --  Theoretical 

During  this  period  a  major  effort  has  been  devoted  to  the  study  of 
hydrogen  kinetics  as  a  tool  for  the  determination  of  electronic  trapping 
levels.  This  has  culminated  in  a  quite  thorough  understanding  of  the 
system  and  its  availability  as  a  tool  for  the  trapping  level  studies.  The 
results  of  the  work  are  reported  in  preprints  1  and  2  attached.  Preprint  1 
deals  with  the  theoretical  aspects  of  the  problem  and  preprint  2  discusses 
the  experimental  results. 

C.  Hydrogen  Kinetics  --  Experimental 


See  preprint  2. 


D.  Devices 

The  problem  of  fabricating  a  reliable  MOSFET  structure  on  rutile  has 
boiled  down  to  three  parts: 

1.  Finding  appropriate  orientations  of  the  crystal  such  that  the 
rutile  substructure  will  not  present  an  inherently  large  surface  charge, 
and  thus  impede  or  eliminate  any  transistor  action  (i.e.,  gain,  or  even 
source-to-drain  conductance).  Transistors  have  been  constructed  successfully 
on  the  (111)-  and  (101)-  surfaces;  however,  devices  constructed  on  (001)- 
and  (110)-  surfaces  have  repeatedly  failed.  Until  our  fabrication  procedures 
have  been  refined  somewhat,  we  cannot  say  whether  this  behavior  is  an 
inherent  property  of  the  surface,  or  a  difference  in  adsorbed  impurities. 

2.  Finding  electrode  material  which  will  not  diffuse  significantly 
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into  the  rutile  while  the  sample  is  heated  to  above  500°C  for  several 
minutes  to  deposit  the  Si02  layer.  In  the  early  transistors,  aluminum 
was  the  only  material  used  for  electrodes.  This  was  marginally  acceptable 
so  long  as  the  source-to-drain  distance  was  ~  200y.  However,  with  the 
recent  use  of  masks  allowing  a  source-to-drain  distance  of  25p,  aluminum 
electrodes  were  found  to  diffuse  quickly  into  the  rutile  and  produce  a  low 
resistance  path  from  source-to-drain.  Other  materials  were  tried,  including 
molybdenum,  copper,  and  gold.  Molybdenum  will  not  diffuse  into  the  rutile 
even  at  much  higher  temperature,  and  thus  appeared  to  be  an  excellent 
prospect.  However,  it  oxidizes  fairly  rapidly  while  the  Si02  is  being  put 
down;  and  worse  yet,  it  appears  to  have  rather  high  mobility  on  the  surface 
of  the  rutile  --  thus,  shorting  out  the  source  and  drain  electrodes.  Copper 
also  oxidizes  rapidly,  forming  layers  of  unstable  and  erratic  conductivity. 
Gold  will  neither  oxidize  nor  diffuse  at  these  temperatures,  but  unfortunately 
it  adheres  to  the  rutile  surface  very  poorly.  Thus,  in  an  attempt  to  make 
use  of  gold  electrodes,  the  following  procedure  was  devised.  Initially, 
semi-transparent  aluminum  electrodes  (about  50-100&)  were  evaporated  onto 
the  rutile,  which  were  then  overcoated  with  evaporated  gold.  This  procedure 
resulted  in  a  strong,  durable  electrode;  A1  diffusion  into  the  Ti02  appeared 
to  be  minimal.  Since  it  was  necessary  to  move  the  sample  when  changing 
from  the  aluminum  evaporating  filament  to  the  gold  evaporating  filament,  a 
highly  stable  masking  device  was  built  which  clamped  the  rutile  sample  over 
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a  1-mil  tungsten  wire.  Thus,  we  now  have  a  highly  reliable,  reproducible 
way  of  producing  a  1-mil  source-to-drain  separation  along  with  rather  good 
electrodes.  Testing  is  proceeding  on  these  devices. 

3.  Finding  a  way  to  clean  the  surface  of  the  rutile  to  eliminate 
mobile  charges  in  the  Si02  and  at  the  Ti02  -  Si02  interface.  A  single 
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monolayer  of  contaminants  can  represent  a  significant  fraction  of  1%  con¬ 
centration  if  it  diffuses  into  a  glass  layer  a  few  thousand  angstroms  thick, 
significant  instabilities  may  develop  and,  in  fact,  this  appears  to  be  a 
major  problem.  For  the  present,  we  are  attempting  to  solve  this  problem 
by  a  combination  of  aqueous  chemical  treatments  and  heat  treating  cycles. 

It  seems  likely,  however,  that  the  final  solution  will  require  fabrication 
of  the  devices  in  the  RF  sputter-etching  system,  which  should  permit  us 
to  fabricate  very  clean  devices.  Installation  of  this  equipment  has  pro¬ 
ceeded  rather  slowly,  since  we  felt  it  was  necessary  to  provide  a  very 
clean  pumping  system,  but  it  is  expected  to  be  operational  within  the  next 

month.  This  should  coincide  well  with  the  availability  of  our  first  locally 
grown  crystals. 


E.  Review  Paper 


1.  Progress  on  the  review  paper  has  been  slow  due  to  the  large  number 
and  complexity  of  published  experimental  results.  However,  a  first  draft 
is  complete  except  for  the  transport  and  sample  growth  sections  and  a 
summary  of  defect  effects.  The  work  involved  in  producing  the  review  has 
given  rewards  in  the  form  of  suggestions  for  useful  experiments  as  well  as 
new  insights  concerning  rutile.  For  example,  we  have  been  able  to  correlate 
optical  and  phonon  data  which,  with  the  optical  selection  rule  tables 
generated  for  the  review,  have  allowed  us  to  conclude  that  the  minimum 
band  gap  in  rutile  is  an  indirect  gap  of  3.02eV,  probably  involving  transitions 
between  the  T  and  M  points  of  the  Brillouin  zone.  These  results  have  been 
submitted  for  publication.  (See  Preprint  III).  These  results  have  oeen 
received  with  considerable  interest  by  another  group  here  at  Utah  which  is 
presently  involved  with  band  structure  calculations  for  rutile.  The  reason 
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for  the  interest  is  that  the  calculated  symmetry  assignments  are  the  only 
experimental  data  available  which  may  be  compared  directly  with  the  band 
structure  calculations. 

F.  Far  Infrared  Spectroscopy 


During  this  period  construction  of  the  spectrometer  has  been  essentially 
completed  with  initial  operation  planned  tor  September.  We  are  setting  up 
the  following  series  of  experiments  for  the  initial  runs. 

The  first  experiment  will  be  to  measure  over  a  wide  range  the  trans¬ 
mission  and  reflectivity  of  the  purest  samples  we  have,  as  provided  by  our 
crystal  growth  effort.  The  purpose  of  these  measurements  is  to  establish 
accurate  values  of  the  transmission  and  reflection  due  to  the  lattice  alone. 
This  data  will  be  stored  on  magnetic  tape. 

The  second  experiment  will  be  to  measure  both  reflection  and  transmission 
of  some  nominally  doped  samples.  This  data  will  be  compared  with  the  pure 
data  to  locate  impurity  induced  absorption  and  reflection.  Once  we  have 
established  the  locations  of  the  impurity  induced  absorption  and  reflection, 
we  will  make  high  resolution  and  high  signal -to-noise  measurements  to 
accurately  determine  the  structure  of  the  impurity  induced  absorption.  Then 
we  will  diffuse  in  impurities  to  determine  whether  the  absorption  is  due 
to  a  given  impurity.  Studies  of  these  absorptions  as  a  function  of  Fermi 
level  (controlled  as  described  in  the  Sections  B  and  C)  will  follow.  Study 
of  absorptions  associated  with  H  and  D  also  have  a  high  priority;  among 
other  things,  such  studies  might  yield  information  on  such  interesting 
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as  the  exact  site  in  the  lattice  of  the  H+. 

Another  experiment  in  this  first  series  will  be  to  look  at  samples 
with  conduction  electrons.  Small  polaron  theory  predicts  that  in  a  perfect 
crystal  there  will  be  no  free  carrier  absorption  below  the  first  interacting 
phonon  energy.  However,  in  imperfect  crystals,  detect  scattering  should 
permit  us  to  measure  free  carrier  absorption. 

fi.  Magnetic  Resonance 

As  mentioned  in  the  previous  (Nov.  1971)  progress  report,  "Vacuum 
Reduction"  does  not  occur  below  1000°C.  However,  some  reduction  does  take 
place  at  higher  temperatures  (~  1200°C).  We  have  analyzed  the  ESR  and 
optical  absorption  spectra  of  rutile  reduced  in  this  fashion.  We  find  an 
hSK  spectrum  similar  to  that  found  by  Chester1  for  vacuum  reduced  material. 
Reoxidation  appears  to  proceed  in  2  stages.  In  the  first  stage,  the  re¬ 
oxidation  proceeds  >  90%  of  the  way  to  completion  as  indicated  by  the  optical 
absorption  spectrum.  However,  the  ESR  spectrum  does  not  change.  Further 
oxidation  is  difficult  but  reduces  the  observed  ESR  spectrum  as  it  proceeds. 

We  find  no  evidence  of  Ti3+  interstitials  and  hence  have  looked  elsewhere 
for  the  source  of  the  spectra.  All  the  ESR  work  was  done  at  4°K  as  the  sample 
is  very  lossy  above  He  temperature.  See  Preprint  IV 

In  an  effort  to  see  if  "vacuum  reduction"  is  due  to  carbon  interstitials 

arising  from  hydrocarbon  contamination  of  dirty  vacuums,  a  sample  was  reduced 
.13 

in  C  0.  Light  optical  coloration  of  the  sample  was  obtained,  but  a  very 
strong  4°K  signal  was  seen  with  g  ,(  =  1.949  and  g  ^  =  1.976.  No  resolved 

hypcrfine  structure  due  to  C  was  seen,  but  there  is  a  suggestion  of  an 
unresolved  splitting.  ENDOk  studies  of  this  possible  splitting  are  planned. 
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As  mentioned  in  Sections  B  and  L  the  Fermi-level  can  be  moved  abouf 
in  the  gap  resulting  in  various  charge  states  for  defects  and  impurities  in 
rutile.  These  slates  have  been  identified  by  ESR  in  the  case  of  Mo.  The 
states  are  Mo  ,  Mo  ,  Mo  ,  and  Mo  with  5+  and  3+  yielding  ESR  spectra. 

This  is  the  first  report  of  Mo3+  in  rutile.  The  spectra  for  Mo5+  have  all 
been  taken  and  are  being  analyzed.  A  paper  on  the  technique  of  Fermi  level 
variation  to  obtain  various  ESR  charge  states  will  be  presented  at  the 
Colloque  Ampere  in  August,  1972. 

H.  Free  Carrier  Optical  Absorption 

A  major  unanswered  question  in  rutile  has  been  the  origin  of  the 
extremely  broad  optical  absorption  peak  centered  at  1.5  y.  It  is  known 
to  be  associated  with  conduction  electrons  since  its  intensity  scales  with 
the  conductivity,  but  the  reasons  for  its  shape  are  much  in  doubt.  The 
most  promising  explanation  for  the  peak  has  been  put  forward  by  Bogomolov,  et  al.2 
who  ascribe  it  to  a  small  polaron  exhibiting  hopping  type  conduction.  A 
key  prediction  of  their  theory  is  that  the  peak  should  shift  to  higher  energies 
as  the  temperature  is  lowered.  They  examined  the  temperature  dependence 
of  the  peak  in  rather  poorly  characterized  samples  and  did  indeed  observe 
a  shift  to  higher  energies.  We  have  now  examined  the  temperature  dependence 
on  well-characterized  samples  and  find  that  for  rutile  with  a  Ti4+  interstitial 
concentration  which  is  relatively  high  compared  to  the  H  concentration, 
there  is  indeed  a  noticeable  shift  of  the  carrier  peak  to  higher  energies 
as  the  temperature  is  lowered  below  77° K.  However,  for  rutile  doped  with 
hydrogen  with  a  nominally  small  Ti4  interstitial  concentration,  no  such 
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peak  shift  is  observed.  Hence,  while  these  experiments  seem  to  confirm 
the  localized  electron  transfer  model  for  the  absorption,  it  is  not  due 
to  an  intrinsic  property  of  the  lattice  but  rather  to  Ti  interstitials. 

In  this  case,  the  peak  shift  corresponds  to  the  trapping  of  electrons  on 

,4+ 

Ti  interstitials.  However,  we  have  also  found  that,  even  at  room  temper¬ 
ature,  rutile  can  be  put  in  a  state  such  that  there  is  a  considerable 
component  of  the  higher  energy  (i.e.,  1.1  eV)  absorption.  This  occurs  in 
crystals  which  have  been  treated  specifically  to  eliminate  all  donors 
(particularly  H)  except  for  Ti  interstitials,  and  in  which  the  Ti  inter¬ 
stitial  concentration  is  fairly  low  --  presumably  of  the  order  of  1/3  the 
density  of  electron  traps.  If  the  absorption  at  1.1  eVis  indeed  due  to 
an  electron  trapped  on  a  Ti  interstitial,  it  is  easy  to  understand  why 
it  has  not  previously  been  seen  at  room  temperature,  since  addition  of 
more  Ti  interstitials  would  result  in  a  strong  free  carrier  absorption, 
which  would  quickly  mask  the  higher  energy  absorption.  Thus,  the  1.1  eV 
peak  can  only  be  observed  over  a  very  narrow  range  of  Ti  interstitial 
concentration,  and  could  easily  be  missed  except  in  a  careful,  systematic 
study  of  absorption  as  a  function  of  Ti  concentration. 

Unfortunately,  the  first  observation  can  be  consistant  with  the 
proposed  model  only  if  the  number  of  electron  states  available  within 
~  -025  eV  (~  kT)  of  the  bottom  of  the  conduction  band  is  comparable  to 
the  concentration  of  Ti  interstitials,  which  seems  rather  improbable. 

Thus,  the  situation  may  be  somewhat  more  complicated.  Work  is  continuing 
on  this  problem;  in  particular,  we  plan  to  do  correlated  studies  of  optical 
absorption  and  AC  conductivity  as  a  function  of  temperature  in  such  crystals. 
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I.  Hydrogen  Diffusion 

As  discussed  in  Section  C,  the  work  on  H  solubility  lias  also  provided 
an  excellent  technique  for  measurement  of  H  diffusion.  The  general  approach 
of  isotope  substitution  probably  provides  the  only  sound  technique  for  such 
measurements  in  systems  where  the  chemical  potential  of  a  given  ion  is 
concentration  dependent.  Although  we  do  not  have  sufficient  data  as  yet 
to  obtain  an  accurate  value  of  the  activation  energy  for  H  diffusion,  we 
have  obtained  some  preliminary  measurements  on  both  doped  and  undoped 
crystals,  and  several  interesting  conclusions  can  already  be  drawn.  For 
example,  we  find  11-diffusion  is  strongly  inhibited  by  the  presence  of  A1 
impurities,  out  much  less  effected  by  Fe.  At  800°C,  the  diffusion  coefficient 
for  an  A1  doped  crystal  is  at  least  10  times  smaller  than  for  a  crystal 
with  comparable  Fe  concentration.  (This  is  in  qualitative  agreement  with 
the  work  of  Hill  ,  although  other  aspects  of  his  work  are  probably  not 
correct).  It  is  tempting  to  try  to  account  for  the  difference  in  diffusion 
coefficients  on  the  basis  of  proton  trapping  near  A1  in  the  lattice,  since 
it  is  known  from  IR  measurements  that  there  is  some  trapping  of  H+  near  A1 
but  apparently  not  in  the  case  of  Fe.  Preliminary  calculations  indicate, 
however,  that  the  trapping  due  to  A1  should  only  produce  a  very  minor  effect. 
Further  study  of  this  problem  is  planned. 

Calculations  indicate  that  the  diffusion  coefficients  for  H+  and  D+ 
at  800°C  should  differ  by  approximately  a  factor  of  2  (D+  being  slower), 
due  to  the  differences  in  zero  point  energy  and  in  the  vibrational  fre¬ 
quencies.  Of.**  measurements,  as  described  in  Section  C,  however,  reveal 
no  measurable  difference.  Indeed,  we  find  that  the  total  concentration 
of  H+  plus  D+  remains  constant  within  experimental  error,  over  the  whole 
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range  of  isotope  substitution,  contrary  to  what  would  be  expected  on  the 
basis  of  the  expected  difference  in  diffusion  coefficients.  Thus  the 
results  of  our  experiments  clearly  illustrate  the  overriding  importance 
of  factors  in  the  driving  force  for  diffusion  other  than  concentration 

gradients  --  specifically,  variation  in  the  electrostatic  potential,  due 
to  space  charges. 

In  summary,  we  believe  we  now  have  access  to  a  technique  which 
permits  meaningful  diffusion  measurements  to  be  made  for  the  first  time 
m  this  material  and  work  is  proceeding  to  exploit  this  technique. 

J.  Trapping  Levels 

As  mentioned  in  Sections  B  and  C,  we  have  developed  a  technique  for 
determining  the  electron  trapping  levels  associated  with  a  given  impurity 
by  means  of  the  dependence  of  the  Fermi  level  on  equilibration  temperature 
in  an  H20  ^ atmosphere.  We  are  currently  applying  this  technique  to  Mo 
doped  samples.  As  noted  in  Section  G,  we  have  identified  the  charge  states 
associated  with  the  optical  absorptions.  Thus,  by  monitoring  the  OH  and 
optical  absorptions  and  using  the  theory  developed  in  Section  B,  w.»  can 
determine  the  thermal  trapping  energies  of  each  state.  As  seen  is  Section 
B.  this  requires  measurement  of  E<H*>,  the  proton  energy.  We  are  currently 
measuring  this  by  the  technique  of  Section  B.  These  same  measurements 
will  also  yield  the  conduction  band  parameters  (the  energy  of  the  conduction 
band  minimum,  Ec,  and  the  effective  density  of  states),  thus  permitting 
placement  of  all  other  measured  energies  relative  to  Ec>  This  measurement 
and  the  Mo  work  is  currently  in  progress  and  results  are  expected  shortly. 


-13- 


K 


17 


We  are  also  preparing  samples  for  similar  work  on  W,  Nb,  V.  and  Fe 
doped  material.  In  fact,  it  seems  that  the  levels  due  to  virtually  any 
impurity  and  charge  state  can  be  mapped  out  in  this  fashion. 
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DEFECT  AND  IMPURITY  THERMODYNAMICS 
IN  RUTILE-LIKE  SYSTEMS* 


John  W.  DcFord  and  0.  W.  Johnson 
Department  of  Physics,  University  of  Utah 
Salt  Lake  City,  Utah  34112 


An  analysis  is  presented  of  the  thermodynamics  of  lattice 
defects  and  H  or  D  in  rutile.  It  is  shown  that  the  solu¬ 
bility  of  these  imperfections  depends  strongly  on  electron 
Fermi  level,  which  in  turn  depends  on  the  concentration  of 
the  imperfections.  A  method  is  presented  for  using  hydrogen 
solubility  to  determine  conduction  band  parameters  and  the 
energies  and  densities  of  electron  traps  in  the  crystal.  The 
conditions  for  obtaining  large  or  small  concentrations  of  H,D 
or  lattice  defects  arc  discussed.  The  use  of  the  method 
presented  depends  on  an  accurate  calibration  of  the  oscillator 
strengths  of  the  H-D  IR  absorption  in  rutile.  Such  a  cali¬ 
bration  is  presented  in  the  following  paper. 


This  research  supported  by  the  Advanced  Research  Projects  Agency 
under  Contract  No.  DAHC-15-71-G8. 


I.  INTRODUCTION 


This  paper  will  discuss  a  theoretical  investigation  of  the  defect 
thermodynamics  of  rutile  and  similar  systems.  Special  attention  will  be 
paid  to  hydrogen  solubility  and  its  application  to  the  determination  of 
electron  trapping  levels  associated  with  impurities  and  lattice  defects. 

The  following  paper1  describes  corresponding  experimental  investigations. 

The  solubility  of  H  or  of  a  lattice  defect  is  complicated  by  the 

fact  that  the  Fermi  level  ■  V*e  -  can  be  widely  varied  (in  principle,  over 

2 

the  entire  band  gap  of  about  3  eV  ).  This  means  that  there  will  not  be 
a  unique  "binding  energy"  for  an  electrically  active  impurity  or  defect 
in  rutile,  since  the  effective  binding  energy  includes  both  the  ion  energy 
(presumably  a  constant)  and  that  of  its  valence  electrons  (which  enter 
the  lattice  at  p  ,  which  is  variable).  Hence,  there  is  no  unique  solubility 
for  such  impurities  or  defects  in  rutile.  We  will  develop  a  complete 
thermodynamic  description  of  the  system,  including  arbitrary  impurities 
and  defects,  subject  to  assumption  that  the  concentrations  are  low  enough 
that  the  perfect  crystal  band  structure  is  not  altered  (an  assumption 
nearly  always  satisfied  in  practice).  We  also  identify  the  parameters 
necessary  to  characterize  the  systems  and  indicate  how  they  may  be  evaluated 
experimentally. 

In  particular  we  are  concerned  with  predicting  and  the  concentrations 
of  H  and  Ti  interstitials  which  will  result  from  various  heat  treatments 
and  atmospheres.  These  two  entities  were  chosen  for  detailed  analysis 
because  H  makes  a  very  convenient,  reversible  dopant  for  probing  electronic 
properties,  and  Ti  interstitials  appear  to  be  the  dominant  defect  and 
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hen^e  will  generally  be  important  in  establishing  equilibrium.  We  will 
also  show  that  direct  and  relatively  simple  measurements  will  suffice  to 
determine  the  trapping  levels  and  concentration  of  any  particular  impurity 
or  lattice  defect.  We  will  derive  expressions  for  the  dependence  of  electron 
Fermi  level  on  these  parameters  and  hence  the  theory  will  allow  prediction 
of  the  treatment  required  to  produce  a  crystal  with  any  desired  pe.  We 
will  show  that  the  key  to  this  technique  lies  in  an  accurate  measurement  of 
H  concentration  and  that  this  can  be  done  by  measurement  of  the  IR  absorption5 
associated  with  H,  once  an  accurate  calibration  factor  is  determined.  The 
following  paper*  includes  such  a  calibration. 

In  spite  of  the  rather  extensive  effort  that  has  been  devoted  to  the 
behavior  of  lattice  defects  and  II  in  the  rutile  system,4,5  there  remains 
considerable  ambiguity  and  much  misunderstanding  as  to  the  conditions  under 
which  hydrogen  and/or  lattice  defects  will  diffuse  into  the  rutile  lattice. 

To  indicate  the  complexity  of  the  system,  we  will  show  that  heating  in  an 
atmosphere  of  1^0  and  0,,  results  in  a  high  concentration  of  H  in  the  lattice, 
while  heating  in  a  dry  1^  atmosphere  generally  results  in  virtually  no  H 
in  the  crystal.  H  is  an  electron  donor  in  Ti02,  and  hence  as  mentioned 
above  the  concept  of  a  unique  "solubility"  for  II  in  rutile  is  meaningless 
until  ue  is  specified.  The  problem  is  further  complicated  by  the  fact  that 
Pe  obviously  depends  on  the  H  concentration  as  well  as  on  the  concentration 
of  all  other  electrically  active  impurities,  some  of  which  are  also  mobile 
at  the  temperatures  involved  in  II  diffusion.  It  is  worth  noting  that  the 
large  band  gap  in  rutile  allows  large  changes  in  pe,  which  can  change 
solubilities  (not  only  of  II,  but  other  species  as  well)  by  many  orders  of 
magnitude.  This  fact  has  been  overlooked  by  many  investigators  and  has 
frequently  resulted  in  misinterpretation  of  results.  It  is,  however, 
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this  dependence  of  H  solubility  on  pe,  and  hence  on  electron  trapping 
levels,  that  makes  it  a  powerful  tool  in  the  analysis  of  impurities,  as 
will  be  shown  below. 

In  Section  II  we  discuss  in  detail  the  theory  of  H  solubility,  the 
determination  of  impurity  trapping  levels  and  conduction  band  density  of 
states  effective  mass,  and  the  distinction  between  the  energies  measured 
by  optical  absorption  and  thermal  energies,  due  to  the  polaron  nature  of 
rutile.**  Possible  applications  of  this  technique  to  various  rutile-impurity 
systems  are  discussed  in  Section  III. 

II.  THEORY 


He  will  now  consider  the  equilibrium  concentration  of  H  in  rutile 
in  the  presence  of  an  atmosphere  of  02,  H2,  and  H20  at  elevated  temper¬ 
atures.  It  is  known  that  rutile  becomes  oxygen  deficient  when  heated  in 

7 

the  presence  of  H2<  This  could  result,  in  principle,  from  the  production 
of  either  0  vacancies  or  Ti  interstitials.  The  preponderance  of  evidence 
supports  the  latter  interpretation,  at  least  at  temperatures  below  1000*C8, 
and  we  will  assume  that  Ti  interstitials  are  the  dominant  defect.  However, 
the  analysis  is  readily  modified  to  allow  for  0  vacancies,  without  quali¬ 
tative  change  in  our  essential  conclusions. 

He  suppose  that  the  system  consists  of  the  following  species: 

HjO,  H2,  02,  Ti02,  <Ti^+>,  <H+>,  e-  where  <Ti*+>  denotes  a  titanium  inter¬ 
stitial  e-  denotes  an  electron,  and  <H+>  denotes  an  interstitial  proton. 

The  valence  of  the  Ti  interstitial  has  received  considerable  attention 
and  it  may  be  that  the  last  word  has  not  yet  been  said  on  this  matter. 
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The  preponderance  of  evidence  is,  however,  that  the  charge  is  "normally" 

4+,  although  there  is  apparently  a  shallow  electron  trap  associated  with 
8 

this  defect.  We  assume  that  this  trap  is  sufficiently  shallow  that  it 
is  essentially  unpopulated  at  the  temperatures  in  question.  The  minor 
modification  of  our  treatment  which  would  be  required  if  this  assumption 
proves  to  be  incorrect  is  again  straightforward. 


The  following  reactions  are  then  possible: 


2H2(g)  ♦  Ti02  ♦+  2H20(g)  ♦  <Ti4+>  ♦  4e' 

(1) 

2H20(g)  4<H+>  ♦  4e~  ♦  02(g) 

(2) 

2H20(g)  <-  2H2(g)  ♦  02(g) 

(3) 

H2(g)  2<H+>  ♦  2e~ 

(4) 

Ti02  «-*  02(g)  ♦  <Ti4+>  ♦  4e" 

(5) 

Eq.  (5)  can  be  obtained  from  Eqs.  (1)  and  (3),  so  these  are  obviously 

not  all  independent,  but  it  is  convenient  to  list  all  possibilities. 

The  condition  for  equilibrium  at  constant  temperature  and  pressure 

is  that  the  total  Gibbs  Free  Energy  for  the  system  be  a  minimum.  We  will 
+  4  + 

assume  that  the  <H  >  and  <Ti  >  do  not  interact ,  except  through  their 

influence  on  pe.  This  is  probably  not  correct  under  all  circumstances, 

4  + 

at  least  for  the  <Ti  >,  since  there  is  some  evidence  for  the  formation 

9 

of  complexes.  Since  there  is  no  detailed  information  available  on  these 
complexes,  it  is  futile  to  attempt  to  treat  them  at  this  point.  Our  treat¬ 
ment  should  be  correct  at  low  concentrations,  however,  and  can  be  readily 
modified  when  experimental  data  on  these  complexes  become  available.  With 
this  assumption,  the  condition  on  the  Gibbs  Free  Energy  can  be  replaced 
by  the  condition  that  the  Gibbs  Free  Energy  for  each  of  the  reactions  (1)  - 
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(5)  be  simultaneously  a  minimum.  Then, 


7 

Z 

i=l 


H.5  dNi  s 

T,P,N 


7 

Z 

i=l 


p.dN.  =  0  , 

i  i  ’ 


(6) 


where  the  derivative  is  taken  with  temperature,  pressure,  and  all  N 
except  constant.  denotes  the  number  of  particles  of  species  i, 
and  vk  is  its  chemical  potential.  (Note  that  throughout  we  have  treated 
the  system  in  terms  of  total  numbers  of  each  specie,  rather  than  concen¬ 
trations.)  The  sum  runs  over  the  7  species  listed  above.  It  is  well 
known  that  p.  can  also  be  written  as 


Vi  = 


3N.  > 


T,V,N 


and  that: 


(7) 


F  =  -kT  log  Z  (g) 

where  Z  is  the  partition  function: 

Z  =  Z  exp[-e./kT],  (9) 

j  J 

where  is  the  energy  of  the  j  n  state  (see  appendix  A),~k  is  Boltzmann's 
constant,  T  is  the  absolute  temperature,  and  the  sum  is  over  all  states 
of  the  system.  Imposing  the  requirement  that  the  total  number  of  each 
chemical  specie  is  constant  leads  to  the  set  of  equations: 

2MH2  +  ^Ti02  "  2mH20  "  U<Ti4+>  ‘  4ye  =  0  (10) 


4p 


<H+> 


(ID 
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2yH20  "  2yH2  ‘  y02  =  °  (12) 

Mm  -  2m  +-2m  =0  (13) 

H2  <H  >  e 


+  u  .  +  4u  -  ii,„.  L 

M  4+  MTi  f) 

<Ti  >  6  1  U2 


=  0 


(14) 


Equations  (10)  -+  (14),  when  accompanied  by  expressions  for  the 

various  chemical  potentials,  completely  determine  the  equilibrium  of  the 

+  4  + 

system.  In  particular,  they  permit  calculation  of  the  <H  >  and  <Ti  > 
concentrations  in  equilibrium  with  a  given  atmosphere  at  any  temperature 
and  pressure. 

The  chemical  potentials  are  either  known  or  readily  calculated  except 
for  ye  and  .  Details  are  given  in  Appendix  A.  The  results  are: 


uH20  -  kT  log  P 
pf|  ■  kT  log  PR 

uo2  *  kT  log  P0 

’  '  E<H*> 


h20  -  \o 


+  kT  log  N 


<H+> 


(15) 

(16) 

(17) 

(18) 


<Ti 


4+. 


=  -  E 


<Ti4+> 


+  kT  log  N 


<Ti4+> 


(19) 


yTi02  =  '  ETi02  +  YTi02 


(20) 


where  the  y’s  are  functions  only  of  temperature  and  are  defined  in  Appendix 

A.  The  P's  are  partial  pressures  in  dynes/cm2.  ue  is  the  electron  Fermi 

level  in  the  crystal  and  depends  on  the  impurities  and  defects  present  (see 

below  and  Appendix  B) .  E  ^nd  ET.n  are  also  defined  in 

<Ti  >  2 
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Appendix  A. 

It  is  convenient  now  to  consider  two  special  cases  in  which 

V  »  PH20  •  P02'  an<1  in  “hich  PH  c  '  P02  ~  PH  ■ 


From  Equations  (10)  +  (20)  we  find:  For  P  »  P  ,  P 

Hz,  0,  1^0 


n<h+>  =  fh2  exp  [ 


H  ,2E<H*>  ■  2Ue  -  ^ 


2kT 


•] 


(21) 


N  +  0^ 

<H  >  2 

-  =  -(p — jr/z  F(T)  exp  (3pe/T;T) 

<Ti4+>  H2 


(22) 


“  2YH  0  "  YTi0  +  3/2  \  ~  E  4+>+  ETi0,+  E<H+> 
where  F(T)  =  exp  [ - - - “ _ 2  <Tl  >  _ 2 

kT 


]•  (23) 


For  P 


h2o 


P°2  "  S 


H2°  °2 
1  <h  +>  =  T~  exp  f — 

P  * 


2YH20  +  4E<H+>  "  4ye 
4kT 


■] 


(24) 


N<H+>  3/4  1/2 

N  =  (p0  ^  (PH  05  g(T)  exp  (3Me/kT),  (25) 

<Ti4+>  2  2 


where 


"  Y 


g  (T)  =  exp  [- 


TiO, 


—  y  -lv 

4  r02  2  yH20 


E<H+>  +  ETi0„ 


-  E 


<Ti4+> 


kT 


]  (26) 


From  Equations  (21)+  (26)  one  secs  that  use  of  a  dry  H2  atmosphere 
isuch  as  would  be  obtained  with  use  of  a  cold  trap)  produces  <Ti4+>  rather 
than  <H+>,  whereas  an  atmosphere  of  H20  +02  produces  <H+>  -  a  somewhat 
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surprising  result.  As  discussed  in  the  following  paper,  this  prediction  is 
borne  out  experimentally. 

From  Equation  (24),  it  is  clear  that  a  measurement,  of  N<H+>  for  given  T, 

P„  n,  and  Pn  determines  (E^„+^  -  y  ).  If  we  can  determine  y  -  E  for  the 
H2(r  02  v  <HT>  e  e  c 

same  sample  we  can  then  combine  the  two  equations  to  get  E<H+>  -  Ec*  Then, 
since  Ec  can  be  determined  from  the  known  work  function^0  we  obtain  E<u+>. 

To  get  the  relation  between  y  and  E  we  need  the  form  of  y  ,  (i.e., 

©  C  6 

its  dependence  on  temperature,  electron  concentration,  and  conduction  band 

parameters).  This  can  be  found  by  an  extension  of  the  work  of  Shockley  and 

Last**,  which  requires  detailed  knowledge  of  the  electron  energy  levels  in 

the  crystal,  including  the  conduction  band  states  and  any  trapping  levels 

lying  in  the  gap.  Experimentally,  many  impurities  can  exist  in  rutile  in 

12 

a  variety  of  different  charge  states  .  Hence,  y  must  be  calculated  for 

0 

trapping  levels  due  both  to  ions  of  two  charge  states  (i.e.,  filled  or  unfilled) 

and  of  many  charge  states  (multiple  trapping  levels).  For  concreteness  we 

consider  two  impurities  --  one  with  only  two  charge  states  and  one  with 

several.  We  suppose  that  each  ionic  charge  state  can  also  occur  in  several 

energy  levels  denoted  by  Eg^  where  the  index  s  refers  to  the  charge  and  the 

index  i  to  the  energy  levels  of  a  particular  charge  state.  We  suppose 

also  that  the  conduction  band  can  be  characterized  by  a  density  of  states 
* 

effective  mass,  m  .  Then  if  n  is  the  number  of  non-valence  electrons  we 

0 

have:  (see  Appendix  B  for  details) 
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n 

e 


N 


-6llcr  -6Eoi 
e  i-e 

»*[■ - -ah-. - ] 


Ee 

i 


-3Eii 


<$•- 


m  0  -BE  . 

N2  Z  se0PeS  Ee  S1 
s=o  i 


E  eBV  E  e‘6Esi 


m 

Z 

s=o 


V  ,2m*kT\3/2f  y\y 

2  {  M2  j  j  1+  y  BCE  HT) 
2tt  iic  o  e7e  c 


(27) 


where : 

Nj  =  number  of  ions  with  only  two  charge  states 
N_  =  number  of  ions  with  m  +  l  charge  states 


V  =  volume  of  the  crystal 
h  =  Planck's  constant  /2tt 
Ec  =  energy  of  the  conduction  band  minimum 
For  Equation  (27),  the  energies  are  measured  relative  to  a  state  in 
which  the  electron  is  removed  infinitely  far  from  the  crystal  (i.e.,  the 
"vacuum  level") . 

Lquations  (24)  and  (27)  can  now  be  used  to  determine  m  and 
(E<H+>  "  Ec).  Consider  a  sample  doped  so  that  there  are  many  trapping  levels 
>  1  eV  below  the  conduction  band  but  very  few  near  it.  (A1  or  Fe  doped 
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crystals  without  <Ti  >  apparently  satisfy  this  criterion,  for  example). 
Then  for  fixed  N<^+>,  the  number  of  conduction  electrons  will  not  change 
significantly  over  a  temperature  range  from  300°K  to  900°K.  Hence,  measure¬ 
ments  of  n£  (the  number  of  conduction  electrons)  at  room  temperature  will 
also  yield  the  number  present  at  the  elevated  temperature  at  which  the 
system  was  equilibrated.  From  Equation  (27)  we  have: 


n£  *  VK(T)3/2I 


(28) 


where: 


*  3/2 

V  /-2m  k' 

K  -(-jf) 


(29) 


and 


1  '  o 


(30) 


(or  the  region  of  low  nc  (i.e.,  E£  -  »  kT) ,  I  can  be  approximated  as: 


-«Bc  '  “e’ 


T  r  m 

I  =  T  e 


(31) 


Equations  (28)  and  (31)  yield  the  following  relation  for  u  : 

e 


2n 


Pe  =  Ec  +  kT  log  I 


K(T)3/2^F  V 


]. 


(32) 


This  provides  the  relation  between  Ue  and  which  we  sought.  Equation 
(24)  can  be  written  in  the  form: 


HO 

Ue  -  E<H*>  *  ^02  ’  ^HjO  *  kT  Io®  C-5— )  -  kT  log  N<h+>  (33) 

P°2 
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-  e<h+>  ’  kT  lo8  (n<h*>)  ♦  h(T)  •  (34) 

Equations  (32)  and  (34)  now  yield: 

"cN<H‘>‘  K  expt-  kf  >  t-CT)  (3S) 

where :  ^ 

L(T)  3  “"yV  exp  [^2-]  (36) 

Hence,  measurements  of  the  product  ncN<H+>  at  two  different  T's  is 
sufficient  to  determine  K  and  (E^^  -  Ec).  (See  Appendix  C  for  details). 

The  work  function  of  rutile  is  known1,0 so  E£  is  determined.  Thus 
E<H*>  can  be  determined*  Once  E<H+>  is  determined,  can  be  found  from 
Equation  (24)  for  any  particular  experiment.  In  cases  where  the  fractional 
occupation  of  any  particular  trapping  level  can  be  determined  from  the 
optical  absorption  spectrum,  Equation  (27)  can  be  used  to  determine  the 
energy  of  the  level.  This  method  will  work  for  levels  separated  by  »  kT, 
which  requirement  can  be  satisfied  in  almost  all  cases.  The  same  information 
can  obviously  also  be  obtained  by  carefully  monitoring  <H+>  concentration  as 
a  function  of  equilibrium  temperature  and  pressures  since  N<H+>  depends  on 
Me  (Eqn.  (24))  and  depends  on  the  energy  of  the  trapping  level  (Eqn.  (27)). 
Again  all  that  is  required  is  that  only  one  level  be  changing  population  at 
a  time  (AE  >>  kT) .  Note  that  this  method  is  often  more  useful  since  it 
eliminates  the  need  to  carefully  resolve  the  absorption  peaks. 

It  must  be  noted  that  the  resulting  energy  differences  between  trapping 
levels  and  the  conduction  band  will  not  be  the  same  as  the  optical  absorption 
energies  since  the  former  is  the  energy  difference  between  "relaxed" 
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(Dolaron)  states  and  the  latter  is  the  energy  of  a  transition  between  a 
"relaxed’1  state  and  an  "unrclaxed"  (i.c. ,  rigid  lattice)  state. 

Once  K  and  E<{1*>  are  known,  it  is  possible  to  systematically  determine 
the  trapping  levels  of  virtually  any  impurity  relative  to  Ec.  By  applying 
Equations  (21)  (23)  to  a  case  in  which  the  trapping  levels  have  been 

determined  and  pg  is  close  enough  to  the  conduction  band  that  it  can  be 
calculated  from  a  measurement  of  nc,  E<Ti4+>  can  be  determined. 

III.  APPLICATIONS 


In  this  section  we  discuss  applications  of  the  theory  of  Section  II 
to  three  sample  problems: 

1.  Minimizing  <Ti4*>  concentration. 

2.  Maximizing  <H  >  without  producing  conduction  electrons  or  <Ti4+>. 

3.  Determination  of  the  trapping  levels  of  an  impurity  with  several 
charge  states,  e.g..  Mo. 

The  first  two  problems  arise  in  the  determination  of  <H+>  oscillator 
strength  as  discussed  in  Ref.  (1).  The  third  is  of  central  importance  in 
any  investigation  of  the  electrical  properties  of  the  crystal. 

To  determine  the  conditions  of  minimum  equilibrium  <Ti4*>  concentration, 
we  combine  Equations  (24)  and  (25)  to  obtain 

N<Ti4+>  =  P02  CXpt(Y02  '  ETi02  +  E<t.4+>  +  YTi02)/kTJexPf-4^e/kT]  (37) 

Examination  of  Eqn.  (37)  makes  it  clear  that  a  large  02  pressure  and  a  high 
Fermi  level  is  required  to  eliminate  <Ti4+>.  The  temperature  dependence  of 
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N  is  also  contained  in  Eqn.  (37),  but  is  more  complicated.  In 

<Ti  >  l 

practice  it  turns  out  that  the  main  concern  is  to  use  a  high  enough 
temperature  that  equilibrium  is  reached  ir.  a  reasonably  short  time.  This 
of  course  is  determined  by  diffusion  times.  Since  practical  considerations 
limit  attainable  values  of  ,  we  require  a  incans  of  producing  a  high  p  . 
This  can  be  done  by  introducing  <H*>.  Hence  we  are  led  to  the  question  of 
how  to  maximize  N<^+>  ,  in  the  presence  of  an  atmosphere. 

The  best  conditions  for  introducing  <H*>  can  be  seen  from  Equation  (24). 
It  is  clear  that  we  need  a  large  pressure.  The  0^  pressure  dependence 
is  weaker  than  in  the  <Ti4  >  case.  Hence  the  optimum  conditions  for  mini- 
mizing  <Ti  >  (while  also  maintaining  a  relatively  high  N<J1+>  level)  are 
to  equilibrate  the  sample  in  high  partial  pressures  of  both  02  and  H20.  As 
discussed  in  Reference  (1),  this  technique  does  indeed  produce  a  large 
<H*>  concentration  and  vanishingly  small  <Ti4*>  concentration. 

The  determination  of  trapping  levels  is  a  particularly  important 
application  of  the  theory.  For  concreteness,  consider  a  sample  doped  with 
an  impurity  having  4  charge  states,  each  of  which  has  only  one  energy  level. 
(In  practice  Mo  has  4  charge  states,  as  determined  by  ESR  measurements  to 
be  published  elsewhere,  but  the  charge  states  have  more  than  one  energy 
level.  This  is  a  non-essential  complication  which  will  be  dealt  with  in  a 
subsequent  publication) .  We  suppose  that  the  experimental  conditions  are 
such  tha*  there  are  no  conduction  electrons.  Then  Equation  (27)  becomes: 


n 

e 


*  ;e  -6<E*-E°-2“e>  .  „ 


;  -^vw 

i  ♦  e  ♦  e  ♦  e 


(38) 


-13- 


34 


where  is  the  number  of  Mo  ions  (assumed  constant)  and  E^  is  the  energy 

of  the  ith  charge  state.  Note  that  0.  +  l  differs  from  Ej  by  the  energy 

of  the  added  electron  plus  any  rearrangement  that  takes  place  among  the 

lower  lying  electrons  on  the  ion.  Hence,  roughly  speaking  the  (i  ♦  l)**1 

level  will  be  partially  filled  if  E.  *  *  E.  +  r  Note  that  this  is  a 

different  type  of  energy  level  scheme  than  is  normally  encountered.  (See 

the  discussion  below  relating  E.  to  E  ) .  The  range  in  values  of  u  for 

which  a  level  is  partially  filled  is  -  kT.  Hence  if  the  levels  are  separated 

by  many  kT ,  only  one  level  will  be  partially  filled  at  a  time.  This  is  the 

key  to  the  method.  Consider  the  case  in  which  E,  ♦  u  z  E_.  Then  le-  A 

l  e  2 

E2  will  be  partially  filled,  E  ♦  y  will  be  »  E. ,  and  E,  ♦  p  will  be 
^  o  c  12© 

<<  Ej.  Hence  Equation  (39)  becomes: 


ne  = 


"Mo'1  *  2c 


-8(E2-Hrue) 


] 


1  ♦  e 


(39) 


Now  if  the  experiment  is  performed  under  conditions  which  keep  the  <Ti4+> 

concentration  low  (high  P.  and  moderate  temperatures)  all  of  the  electrons 

2 

must  be  provided  by  <H»>.  Hence  we  can  set  n  *  n ...  .  Rcferri  e  to 

e  ’ 

Equation  (24)  we  sec  thn-t  measurement  of  N^„  .  determines  u  .  In  fact: 

<H*>  e 


Bu, 


<V 

v* 


cxp[- 


r0„  "2yH„0  +  4E<H*> 


2 

ITT 


<H+> 


'V 

^7 


J(T) 


N 


<n+> 


(40) 
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Combining  Equations  (39)  and  (40)  yields: 


S(e2-e1) 


1  1-X  J 


,2NMo  ■  N<H.>,  (V 

1n<H.>-  NMo  1 


(PH205 


J(T) 


V 


J(T) 

N<1.*> 


(41) 


N, 


where  X 


Mo 


of  N, 


N<H+> 
vs  T 


.  Now  X  can  be  determined  at  any  point  from  a  plot 
,<H+>  * eqUi i »  wbere  TeqUil  is  the  e(luilibration  temperature.  Such  a 

plot  is  shown  schematically  in  Fig.  1.  By  measuring  the  height  of  a  step 
we  obtain  N^,  and  hence  we  can  determine  X  at  any  point  on  the  curve. 

Thus  E2  -  Ej  is  determined.  A  similar  treatment  of  the  other  steps  will 
yield  all  the  energy  differences.  (In  practice,  it  may  be  more  feasible 
to  vary  parameters  other  than  Tequil»  but  similar  results  obtain.) 

There  remains  only  the  question  of  where  the  levels  E^  lie  with  respect 
to  the  conduction  band.  The  pertinent  question  is  how  much  energy  is 
required  to  convert  from  E^  to  E2  plus  an  electron  in  the  conduction  band? 
Let  that  energy  be  A.  Then 


E_  ♦  A  =  E_  ♦  E 
3  2  c 


Thus: 


A  =  E_  ♦  E-  E_ 

2  c  3. 


Now  E.  -  E,  is  known  as  shown  above.  E  is  known  as  discussed  in  section  II. 
Hence  A  is  determined  and  all  levels  are  placed  with  respect  to  the  conduction 
band. 
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This  technique  is  currently  being  applied  to  a  variety  of  transition 
metal  impurities  in  Tit^;  results  will  be  reported  in  future  publications. 

IV.  SUMMARY 


We  have  presented  a  method  for  determining  conduction  band  density  of 
states  effective  mass  and  solubilities,  binding  energies,  and  trapping  levels 
of  various  impurities  in  rutile  and  similar  systems.  The  method  relies  on 
an  accurate  determination  of  the  <H+>  oscillator  strength.  Such  a  determination 
is  presented  in  the  following  paper,  Ref.  1. 
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Appendix  A 

The  partition  functions  of  gases  can  be  factored  into  translation, 
rotation,  vibration  and  electronic  parts1 The  translation  part  can  be 
treated  classically  at  the  temperatures  of  interest  and  is 

Ztr«n,  '  jz  t*"")3'2  (1A) 

In  the  cases  of  H2,  02,  f^O,  only  the  ground  electronic  state  contributes. 
Hence: ^ 

-r  AT 

Zelect  "  V  <2A> 

where  rQ  is  the  degeneracy  of  the  ground  state  and  is  the  energy  of 
the  ground  state.  The  vibrational  part  involves  a  sum  over  harmonic  oscillator 
energy  levels  for  each  vibrational  mode.  In  the  cases  of  H 2  and  02,  only 
the  ground  vibrational  state  contributes.  Hence  the  zero  point  energy  and 
the  electronic  ground  state  energy  combine  to  give  the  dissociation  energy 
and  thus: 

VkT 

ZelectZ»lb  ’  V  <»> 

where  ep  is  the  dissociation  energy.  In  the  case  of  H20,  the  first  excited 
states  of  the  3  nor?,al  modes  also  contribute  and  yield:1* 

eD/kT  1 

Zelect  Zvib  '  V  (  -»&>  (  .sis3) 

1-e  T  ! -e  T 


-5404 
1-e  T“ 
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Finally,  the  rotational  parts  are  complicated  by  the  symmetries  of  the 

13 

nuclei.  However,  they  are  well  understood  and  yield  (at  -  300°C  and  above) 


ZD  .  a  2.431  x  10  T 

Roto2 


3.358  x  10"2  T2 


(4A') 


(4A") 


ZRotH  "  2'424  *  10  T 
H2 


Combining  Equations  (1A)  (4A")with  Equations  (7)  and  (8) 
yields : 


*  5.896  x  104  l 

Yn  =  ♦  1.380  x  10_1°  T  log  {3.421  x  10°  e  T  T  1  } 

U2 


P  *  kT  log  P  -  y0 
2  u2  U2 


.  1.111  x  10J  . 

Yh  0  ■  1.380  x  10'10  T  log  (6.645  x  10*  e  T  T4 


^  -2296^  (  -5253^  (  -5404^ 

1-e  T  1-e  T  1-e  T 


mh2o  3  kT  log  ph2o  "  yh2o 


(5A* ) 


(6A') 


,  5.195  x  104  l 

1.380  x  10' 10  T  log  {5.330  x  10^  e  T  1*  } 
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=  kT  log  P 
2  H2 


(7A  * ) 


The  partition  functiore  for  <H+>  and  <Ti4+>  are  found  by 
assuming  that  each  has  but  one  state  in  the  crystal.  This  is 
not  quite  true  in  the  case  of  <H+>  as  snown  by  slightly 
displaced,  impurity-associated  absorption  peaks15.  This 
assumption  is  a  good  approximation,  however,  and  can  easily 
be  removed  if  necessary.  We  chose  as  a  standard  state  --  the 
zero  of  energy  --  the  state  of  separated  H  atoms,  0  atoms,  and 
Ti  atoms.  We  will,  however,  measure  electron  energies  relative 
to  the  state  in  which  they  are  infinitely  separated  from  every¬ 
thing  ("vacuum"  level). 

Consider  <H+>.  Starting  with  the  standard  state  we  produce 
<H  >  as  follows:  First  remove  the  electron.  This  requires 
energy  ER.  Next,  insert  the  electron  in  the  crystal.  This 
takes  energy  u  .  Then  insert  the  proton  in  the  crystal.  This 
requires  energy  Ep.  Thus  the  total  energy  required  to  produce 
<H+>  from  the  standard  state  is: 


Ep  *  ER  .  li, 

(Note  that  both  Ep  and  pe  are  negative  numbers ) . 

We  chose  totreat  the  electrons  separately.  Hence,  let 

(8A) 


Consider  <Ti4+>. 


.4+ 

To  produce  a  <Ti  >  we  must  remove  4 
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electrons .  This  takes  an  energy  E^  +  +  E^  +  E^  We  then 

place  the  4  electrons  in  the  crystal .  This  requires  an  energy 

4+ 

4pg.  We  then  insert  the  Ti  in  the  lattice.  This  takes  an 

4+ 

energy  Ep.  Hence  the  total  energy  required  to  make  a  <Ti  > 
is 

E1  +  E2  +  E3  +  E4  +  4pe  +  Et  > 

Again  we  treat  the  electrons  separately  and  define 


(9A) 

“E<Ti4+>  =  Ei  +  E2  +  E3  +  E4  +  ‘ 

We  can  now  calculate  the  partition  functions  from  Equations 
(7)  -(9).  We  have 


N! 


F  =  -  kT  {N  pjr  -  N  log  N  +  N} 


p  =  -  E  +  kT  log  N 


Hence : 


(10A) 

P  =  -E  +  kT  log  N 
<H+>  <H  >  <H  > 


<Ti 


4+. 


=  -E  .  +  kT  log  N  . 

_.4+  6  „.4+ 

<Ti  >  <Ti  > 


(HA) 
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We  can  calculate  in  the  solid  state  as  follows, 

We  have 


C  =  -T  ) 
v  aT7 ' 


V 


Hence : 


T  C  T  . 

1  T  dT  "  -  1 

o  1  O  ?T  v  3T  V 


dT 


3F  ,  3F  i 

3t5v  I  +  aT^y  I 

V  t  V  q 


§>v  lT  -  s<°> 


where  S(0)  is  the  entropy  at  0°K.  Let 


(12A) 


o 


where  Cy  is  the  specific  heat  at  constant  volume.  Then 


T 

I 

o 


8F.  AT 

W\  dT 


T 

/  [A  -  S(0)]  dT 
o 


(13A) 

T 

F(T,V)  -  F(0,V)  =  -  S(0)T  +  /  AdT  . 

o 

Now  let  -E™...  be  the  energy  of  the  ground  state  of  a  TiO„ 
2  2 

molecule  in  the  crystal  relative  to  the  standard  state. 
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Then, 


( 14  A) 


F(0,V)  =  -  NIL,.- 

2 

where  N  is  the  number  of  TiC>2  molecules  in  the  solid.  Letting 

(15A) 


where  Cy  is  the  specific  heat  per  molecule,  and  noting  that  S(0) 
we  find 


^Tit^  "  "^Ti02  +  YTi02 
where 


T 

-/  dT* 
o 


T' 

/ 


C  (T") 


dT" 


(16A) 


(17A) 
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Appendix  B 


The  chemical  potential  for  N  electrons  is 


-0e. 

Z  (N)  =  I  e  1 
i 


(IB) 


subject  to  the  restriction  that 


(2B) 


In.  =  N 
.  1 

l 


where  the  sum  is  over  all  states  of  the  system.  The  restriction 
can  be  removed  in  the  usual  way1**  by  calculating  instead  the 


function  zg  given  by 


-Y  N' 

zo  =  ZZ.(N')e 
N' 


(3B) 


and  choosing  y  so  that  Z  (N')e 
©  © 

Then 


-Y  N* 


has  a  maximum  at  N'  =  N. 


log  ze  =  log  Ze(N)  -  yeN 


F  =  -kT  log  Z  (N)  =  -kT  log  z  -kT  y  N 
©  6  6 


Thus, 


\i  =  -kT  y 
e  'e 


(4B) 


(SB) 


(6B) 


The  condition  on  yg  leads  readily  to  the  relation 
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* 


N  = 


d  log  z 
*  e 

dYe 


(7B) 


We  will  assume  that  the  impurities  are  independent  (except 
for  their  effect  on  Mg) .  Then  the  partition  function  factors 
into  terms  for  the  conduction  band  and  for  each  of  the  impurities. 
Thus, 

log  ze  =  log  Zj  ♦  log  z2  ♦  log  zcQnd  .  (8B) 

Now,  since  each  impurity  ion  is  assumed  independent,  the  partition 
function  for  a  given  kind  of  ion  factors  into  terms  for  each 
individual  ion.  This  means  that 
N. 

zi  =  W  1  (9B) 

where  z  is  the  partition  function  for  a  single  ion  and  Ni  is  the 
number  of  ions.  Thus, 

m  6m  s  -6E  . 

log  Zj  =  Nj  log  E  e  Ee  SJ  (10B) 

s=0  j 

Equations  (6B) ,  (7B),  (8B),  and  (10B)  can  then  be  combined  to 
yield  Equation  (27).  The  conduction  band  term  in  (27)  arises 
from  ordinary  Fermi-Dirac  statistics^. 
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Appendix  C 


To  determine  the  values  of  K  and  E  ,  imagine  measurements 

<H*> 

of  N  +  n£  at  two  temperatures,  T  and  T.  ♦  A.  Equation  (35) 
<H  ^ 

then  yields 


,L(T2)  tncN<H*>’l. 


P  .  c  _  k(Tl  A)T1  r“v'2'  ~  <H-  , 

E<H  >  c  A  ln  [LCfTT  (n  N  )  1 

1  c  <H+>  2 


(1C) 


(n  N  ) 
v  c  ♦ 

<H  >  1 


K  = 


<H 

"LCfjT 


L(T1)(ncN..>  > 


<H  > 


l(T2) (ncN<„*>,l 


T 

t  1  +  A. 

(  z  ^ 


(2C) 


where  T_  =  T  ♦  A  and  (n  N  )  refers  to  the  value  at  T  , 

C  <H  >  1  1 

and  (n  N  )  refers  to  the  value  at  T_. 

<H  >  1 
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"An  Experimental  Technique  for  the  Precise  Determination 
of  II  and  D  Concentration  in  Rutile  (Ti(\)."* 


0.  W.  Johnson  and  John  DeFord 
Department  of  Physics,  University  of  Utah 
Salt  Lake  City,  Utah  84112 


Results  of  a  precise  determination  of  the  IR  absorption 
strength  of  H  and  D  in  rutile  are  presented,  which  permit 
determination  of  H  or  D  concentration  to  an  accuracy  of  3%. 
Used  in  conjunction  with  the  results  of  the  theoretical 
analysis  of  the  preceding  paper1,  this  technique  makes  pos¬ 
sible  the  accurate  determination  of  the  energy  and  density 
of  all  electron  trapping  levels  in  Ti02  crystals.  Several 
previously  unreported  satellite  spectra  for  both  II  and  I)  arc 
described  and  possible  identifications  arc  discussed.  A 
technique  for  measuring  diffusion  parameters  of  II  and  D 
undcj  conditions  of  strict  chemical  equilibrium,  thus  avoiding 
the  large  departures  from  Fick’s  law  which  usually  plague 
such  measurements  in  materials  such  as  TiO, ,  is  also  described. 
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I.  Introduction 

The  preceding  paper *  outlines  a  method  of  analysis  by  means  of  which 

the  complete  equilibrium  "defect  thermodynamics"  of  u  system  such  as  rutile 

can  be  derived  from  measurements  of  the  solubility  of  H  as  a  function  of 

equilibration  temperature  and  partial  pressures.  By  this  we  mean  that  the 

electron  energy  levels  and  densities  of  all  electron  traps  can  be  determined 

for  a  crystal  containing  unknown  impurities,  or  alternatively,  the  electronic 

structure  and  H  solubility  of  a  crystal  of  known  doping  can  be  calculated. 

In  addition,  the  analysis  permits  calculation  of  the  equilibrium  concentration 

of  Ti  interstitials  (which  we  assume  to  be  the  dominant  lattice  defect  and 

4+ 

will  hereafter  denote  as  <Ti  >}and  an  independent  measurement  of  the  conduction 
band  effective  mass  (or  the  conduction  band  effective  density  of  states). 

The  importance  of  these  parameters  in  characterizing  a  semiconductor  for 
device  applications  (either  electronic  or  electro-optical)  or  in  achieving 
a  fundamental  understanding  of  the  material,  is  sufficiently  obvious  as  to 
require  no  further  comment. 

Utilization  of  this  analysis,  however,  requires  a  sensitive,  accurate 
means  of  monitoring  the  concentration  of  <II+>  (II  interstitials)  in  the  crystal. 
Furthermore ,  the  measurement  technique  must  be  non-destructive,  since  the 
accuracy  required  demands  that  successive  measurements  be  performed  on  the 
same  specimen  to  eliminate  error  due  to  variation  in  sample  purity,  which  is 

unavoidable  with  presently  available  crystals.  Measurement  of  the  IR 

♦  2 
absorption  associated  with  <11  >  satisfies  all  Of  these  requirements. 

However,  utilization  of  this  technique  requires  rather  precise  determination 
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of  the  "absorption  strength"  of  the  <H+>.  Theoretical  calculations  of 
sufficient  accuracy  arc  not  feasible,  because  of  the  large  dielectric 
constant  and  strong  anisotropy  of  rutile,  uncertainties  as  to  the  <ll*> 

7 

site  in  the  lattice,  etc.  A  previous  paper”  gave  an  approximate  value 
for  the  absorption  strength,  based  on  changes  in  mass  of  a  crystal  into 
which  H  had  been  diffused,  but  the  value  listed  was  not  sufficiently 
accurate  for  our  purpose.  In  this  paper  we  present  the  results  of  a 
refinement  of  this  measurement,  yielding  accurate  values  of  the  absorption 
strength  for  both  <H*>  and  <D+>,  together  with  a  detailed  discussion  of 
techniques  developed  to  overcome  the  rather  severe  experimental  difficulties 
inherent  in  such  a  measurement.  Work  is  presently  in  progress  to  apply 
these  results  and  the  analysis  in  the  preceding  paper1  to  a  variety  of 
rutile-impurity  systems. 


II.  Description  of  Experiment 

The  basic  requirement  of  the  experiment  is  very  simple;  to  determine 
precisely  the  concentration  of  <H+>  required  to  produce  a  given  IR  absorption 
Several  techniques  were  considered,  including  various  chemical  determinations 
mass  spectrometer  measurements,  NMR,  and  others.  For  various  reasons,  it 
was  decided  that  direct  measurement  of  the  mass  changes  associated  with 
changes  of  <H*>  concentration  were  most  likely  to  yield  the  required  accuracy 
In  practice,  it  was  found  that  this  straightforward  approach  was  inadequate, 
for  reason*  discussed  below,  and  it  was  necessary  instead  to  measure  the  mass 
change  associated  with  the  replacement  of  <H*>  in  the  crystal  with  <D*> 
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(or  conversely),  in  such  a  way  that  the  total  concentration  of  <H+>  plus 
<D  >  remained  essentially  constant.  Thus,  the  measurements  consisted  of 
weighing  a  sample,  equilibrated  at  appropriate  H,0  partial  pressure  and 
temperature,  measuring  the  IR  absorption  due  to  the  <H+>,  then  heating  the 
crystal  to  the  same  temperature,  wita  the  1I20  replaced  by  0,0,  reweighing 
the  crystal  and  measuring  the  IR  absorption  due  to  both  <D+>  and  <H+>.  By 
cycling  the  crystal  several  times  and  checking  at  various  points  in  the 
cycle,  it  was  established  that  the  total  concentration  of  <H+>  plus  <D*> 
was  indeed  constant  (except  as  noted  below),  and  the  absorption  strengths 
of  both  the  <H+>  and  <D+>  were  readily  obtained.  A  variety  of  complications 
arose  in  the  course  of  the  experiment,  as  discussed  below,  but  the  basic 
experiment  consisted  of  this  simple  sequence  of  measurements. 

III.  Sample  Preparations  and  Heat  Treatment 

As  shown  in  the  preceding  paper1,  the  <H+>  solubility  is  strongly 

dependent  on  the  den-tty  of  acceptor  levels  in  the  crystal,  and  because  of 

the  very  small  ma.^s  changes  expected,  it  was  necessary  to  utilize  crystals 

containing  as  high  a  concentration  of  acceptors  as  possible.  (An  alternative 

approach  would  be  to  H-dope  without  acceptors,  allowing  the  electrons  from 

the  H  donors  to  populate  the  conduction  band;  however,  the  conduction  electron 

absorption  is  much  stronger  than,  and  overlaps  the  <H+>  absorption  peak2, 

making  accurate  measurements  unfeasible.  Fe  and  A1  were  chosen  as  the  most 

appropriate  acceptor  impurities,  since  both  have  fairly  high  solubility  and 

were  known  to  diffuse  quire  readily.  Nominally  pure  crystals  were  doped 

by  diffusion  in  sealed,  evacuated  quartz  ampoules,  using  A1C1-  and  FeCl 

3  3 
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as  the  dopants.  It  was  found  that  most  of  the  dopant  (~  80°u)  diffused  into 
the  rutile,  after  heating  to  -  1000°C  for  10  to  100  hours,  depending  on 
sample  dimensions,  and  that  reasonable  good  control  of  acceptor  concentration 
could  be  obtained  by  careful  weighing  of  the  dopant.  The  solubility  of  Fe 
was  found  to  be  somewhat  higher  than  that  of  Al.  Not  surprisingly,  in  the 
light  of  the  discussions  in  Ref.  1,  both  were  found  to  depend  strongly  on 
Me  (and  hence,  on  the  atmosphere  and  temperature  of  heat  treatment).  In  an 
atmosphere  consisting  of  22  Torr  H20  and  100  Torr  02,  at  ~  800°C,  the  Fe 
solubility  was  estimated  to  be  -  6  x  1019/cm3,  and  Al  solubility  was  somewhat 
smaller.  When  the  H20  was  removed  from  the  atmosphere  (for  example,  by  use 
of  a  liquid  nitrogen  trap),  the  solubilit)  decreased  dramatically  in  both 
cases,  resulting  in  diffusion  of  the  Fe  or  Al  to  the  surface  of  the  crystal 
and  formation  of  visible  precipitates  inside  the  crystal.  Since  formation 
of  these  precipitates  involved  a  change  in  the  effective  valence  state  of 
the  dopants  (presumably  both  ended  up  in  +3  states,  while  they  had  previously 
occupied  lattice  sites  with  an  "effective"  valence  of  +4),  this  process  was 
accompanied  by  loss  of  oxygen  from  the  crystal  and  a  corresponding  mass 
decrease,  much  larger  than  the  expected  change  associated  with  <H+>  loss. 

While  this  process  has  not  been  studied  in  detail,  it  is  clear  that  the 
solubility  of  acceptor  impurities  depends  strongly  on  and  that  it  is  not 
possible  to  maintain  a  high  concentration  in  the  absence  of  donor  impurities. 
Thus,  it  was  not  possible  to  perform  the  desired  measurements  with  sufficient 
accuracy  simply  by  diffusion  of  H  in  and  out  as  originally  intended.  To 
circumvent  this  problem,  it  was  necessary  to  perform  the  experiment  in  such 
a  way  as  to  always  maintain  a  constant  (equilibrium)  concentration  of  <H+> 
or  an  entity  chemically  equivalent  to  <H+>,  i.e.,  <D+>. 
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A  further  complication  in  such  an  experiment  is  that  the  concentration 

of  donor  impurities  or  defects  other  tl m  <H+>  or  <D+>  must  be  kept  constant 

and  at  a  very  low  level.  Introduction  of  donor  impurities  can  be  effectively 

eliminated  by  careful  attention  to  possible  sources  of  contamination  and  by 

maintaining  strongly  oxidizing  conditions  during  heat  treatment.  (Hydrocarbons 

from  vacuum  pump  oil,  0-rings,  etc.,  are  a  major  source  of  contamination;  an 

oxidizing  atmosphere  renders  these  innocuous).  All  heat  treatments  were 

performed  in  clean  Vycor  tubes,  in  atmospheres  of  0 2  and  either  Ho0  or  Do0, 

and  we  found  no  evidence  of  contamination,  even  after  as  many  as  20  runs  on 

individual  specimens.  The  problem  of  lattice  defects  is  more  complicated. 

Since  the  introduction  or  removal  of  <Ti  >  (or  0-vacancies)  again  involves 

large  changes  in  mass,  it  was  essential  that  these  defects  either  be  eliminated 

or  at  least  maintained  at  a  constant  level.  After  initial  doping  with  Fe  or 

Al,  the  concentration  of  defects  was  quite  high.  (We  assume  that  the  defects 

are  predominantly  <Ti4+>  or  possible  complexes  of  interstitials,  as  mentioned 

in  Ref.  1.  However,  nothing  in  our  experiment  would  distinguish  between 

these  defects  and  0-vacancies  and  none  of  our  conclusions  would  be  changed 

if  0-vacancies  were  involved.)  This  high  initial  concentration  of  defects 

was  reduced  below  measurable  levels  by  prolonged  treatment  in  an  Hn0  +  0., 

atmosphere  at  temperatures  of  800°C  or  more,  as  discussed  below.  This  treatment 

was  accompanied  by  a  mass  increase  appropriate  to  removal  of  a  concentration 
.  4  + 

of  <Ti  >  corresponding  roughly  to  that  required  to  compensate  the  acceptor 
impurities  which  had  been  introduced;  an  appropriate  simultaneous  increase 
m  <H  >  concentration  was  also  observed  by  IR  absorption. 

Disc  shaped  samples  8  to  20  mm  diameter  and  2  to  10  mm  thick,  oriented 
with  the  flat  surface  i  to  the  c-axis  (to  eliminate  problems  associated 
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With  the  dichroism  of  the  <H+>  or  <D+>  optical  absorption2)  were  used, 
the  dimensions  being  chosen  to  optimize  the  sensitivity  of  the  various 
measurements  and  to  provide  reasonable  diffusion  times  at  the  temperatures 
used.  After  Fe  or  A1  diffusion  and  the  H20  pretreatment  outlined  above,  the 
crystals  were  carefully  polished  on  all  surfaces.  This  step  is  particularly 
important,  since  a  rough  surface  is  very  susceptible  to  chipping  on  a 
microscopic  scale,  and  is  difficult  to  keep  clean.  The  samples  were  checked 
for  homogeneity  by  measuring  the  <H+>  -  <D+>  concentration  in  various  parts 
of  the  sample.  The  A1  doped  crystals  were  found  to  k.  rather  inhomogeneous; 
the  A1  concentration  appeared  to  be  as  much  as  25%  higher  near  the  edges  of 
the  crystal  than  at  the  center.  No  significant  inhomogeneities  were  found 
in  the  Fe  doped  crystals;  <H+>  concentration  was  constant  to  within  2% 
throughout  the  crystal.  For  this  reason,  and  because  of  the  higher  solubility, 
most  of  the  measurements  were  made  on  Fe  doped  material;  results  on  A1 -doped 

crystals  were  consistent  with  those  for  Fe-doped  crystals  within  experi¬ 
mental  error. 


IV.  Optical  Absorption  Measurements 


IR  absorption  measurements  were  made  on  a  Beckman  IR-12.  The  peak 
optical  density  of  the  <H+>  or  <D+>  absorption  line  was  generally  too  high 
for  direct  measurement  (due  to  the  limitations  imposed  by  the  sensitivity 
of  the  mass  measurements  discussed  below),  so  measurements  were  made  at 
vario.-s  fixed  wavelengths  on  the  "tails”  of  the  absorption  lines,  and 
the  total  absorption  strength  (i.e.,  the  area  under  the  absorption  peak) 
was  inferred  from  this  by  careful  comparison  of  absorption  at  these  wave- 
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lengths  with  the  total  absorption  in  thin  slices  of  the  crystal.  It  was 
necessary  to  do  this  separately  for  each  specimen,  since  the  line  shape  and 
background  absorption  varied  somewhat  from  sample  to  sample,  and  was,  of 
course,  quite  different  for  Fe  and  A1  doped  samples.2  Various  checks 
indicated  that  these  determinations  were  reliable  to  -  ±  L'o. 

Since  the  <H  >  and  <D  >  absorption  line  petitions  and  line  widths 
are  quite  sensitive  to  temperature,  it  was  necessary  to  utilize  a  control led - 
temperature  sample  holder.  All  measurements  reported  herein  were  taken  at 
35°C. 


V.  Mass  Measurements 

The  crystals  were  weighed  on  a  modified  Mettler  M-5  microbalance. 

Optimum  sample  size  was  in  the  range  from  5  to  20  grams  (the  maximum 
capacity  of  the  balance)  and  the  corresponding  mass  changes  expected  were 
-  70-300  y-grams  for  crystals  with  maximum  Fe  concentration.  After  con¬ 
siderable  effort,  we  were  able  to  obtain  sensitivity  and  reproducibility 
of  betcer  than  ±  2yg  on  a  5g  sample,  which,  with  adequate  statistics, 
should  permit  an  overall  accuracy  of  1  or  2%,  which  was  the  goal  of  the 
experiment.  The  absolute  accuracy  of  the  mass  determinations  was  substantially 
poorer,  of  course,  but  was  of  no  particular  importance  in  this  experiment. 

To  achieve  this  degree  of  reproducibility,  rather  elaborate  pre¬ 
cautions  were  necessary.  As  previously  mentioned,  careful  polishing  of 
the  ciystals  was  essential  and  extreme  caution  in  handling  was  still  necessary. 
Plastic  tipped  forceps  were  used  to  avoid  chipping  and  scratching,  and 
sliding  the  crystal  against  the  walls  of  the  Vycor  heat-treating  tube  was 
carefu’ly  avoided.  Ordinary  cleaning  techniques  proved  to  be  adequate, 
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provided  all  lint  and  dust  was  subsequently  removed  with  an  antistatic- 
brush.  Adsorption  of  moisture,  etc.,  on  the  surface  of  the  crystal  was 
also  minimized  by  the  careful  polishing  and  was  found  not  to  be  a  problem. 

Appropriate  measures  were  taken  to  isolate  the  microbalance  from 
vibration  and  it  was  located  in  a  room  chosen  to  minimize  temperature 
fluctuations  (the  temperature  was  stable  to  ±  1/2°C,  with  negligible  short 
term  fluctuations).  Nevertheless,  temperature  stability  of  the  balance 
was  still  a  major  problem.  The  balance  was  heavily  insulated  and  the 
front  panel  was  carefully  maintained  at  a  temperature  ~  1°C  above  ambient 
by  heating  elements  cemented  to  the  panel.  The  heating  elements  were  turned 
off  while  the  balance  was  in  use;  this  procedure  compensated  for  the 
temperature  rise  which  otherwise  would  have  resulted  from  the  proximity 
of  the  operator.  It  was  also  found  to  be  necessary  to  provide  forced- 
air  cooling  for  the  light  used  to  illuminate  the  balance  scale.  Remote 
handling  of  the  samples  inside  the  balance  chamber  was  provided.  Intro¬ 
duction  of  samples  into  the  balance  chamber  produced  fluctuations  requiring 
an  equilibration  time  of  ~  1/2  hour.  Appropriate  corrections  to  compensate 
for  variations  in  atmosphere  density  were  necessary,  since  the  density  of 
Ti02  is  significantly  less  than  that  of  the  balance  weights.  With  these 
precautions,  the  measured  masses  of  both  the  rutile  specimens  and  stainless 
steel  standards  were  reproducible  to  ±  3  pg  for  single  readings  and  signifi¬ 
cantly  better  if  several  readings  were  taken. 

VI.  Results  and  Discussion 

A  large  number  of  measurements  were  made  on  a  variety  of  crystals 
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(primarily  Fe-doped,  as  discussed  previously)  and  under  differing  equilibration 

conditions.  The  data  reported  here  are  for  crystals  equilibrated  at  800‘C, 

in  -  300  Torr  02  and  1500  Torr  1^0  or  D^O.  The  exact  0,,  pressure  could  not 

be  determined,  but  was  the  same  from  run  to  run  to  within  SI. .  The  crystals 

were  placed  in  Vycor  capsules  with  ~  0.1  cm3  H20  or  D.,0  which  was  then  frozen. 

Next  the  capsules  were  evacuated,  flushed  and  filled  with  0,  to  a  pressure  of 

100  Torr,  at  which  point  the  capsules  were  sealed.  The  capsule  was  then 

placed  in  a  2-zone  furnace,  with  the  end  containing  the  crystal  maintained 

at  800  C  and  the  opposite  end  at  120°C,  resulting  in  the  partial  pressures 

indicated  above.  At  02  pressures  significantly  below  this  value,  we  observed 

mass  changes  as  a  function  of  02  pressure  which  indicated  a  changing  and  non- 

negligible  concentration  of  <Ti4+>.  Data  taken  in  this  range  and  at  different 

temperatures  can  be  used  to  obtain  the  solubility  parameters  for  <Ti4+>  and 

will  be  reported  in  a  future  publication,  but  for  present  purposes  it  is 

significant  to  note  that  data  obtained  at  02  and  H20  pressures  below  300  and 

1500  Torr,  respectively,  are  consistent  with  the  results  reported  here  only 
.  4+ 

if  <Ti  >  concentration  is  taken  into  account  in  calculating  both  the  mass 

changes  and  the  donor  density.  At  the  pressures  and  temperature  indicated 

above,  our  data  show  that  the  Ti  interstitial  concentration  is  below  detectable 

limits  (2.5  x  1016/cm3).  If  the  <Ti4+>  concentration  is  below  detectable 

limits,  small  unavoidable  changes  in  partial  pressures  and  temperatures 

produce  minor  changes  in  <H  >  and  <D  >  concentration,  which  is  properly 

accounted  for  by  our  analysis,  but  no  detectable  mass  changes  due  to  changing 
.4+ 

<Ti  >  concentration;  hence,  these  data  are  far  more  accurate  than  would  be 
possible  in  the  presence  of  significant  <Ti4+>.  Thus,  only  a  few  runs 
were  required  to  obtain  sufficient  accuracy  after  the  appropriate  operating 
conditions  were  determined.  Since  our  main  concern  here  is  with  the  "absorption 
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strength"  of  <H  >  and  <D  >,  we  report  only  data  for  crystals  with  negligible 
4* 

<Ti  >.  All  of  our  data  (>  100  runs),  however,  are  consistent  with  that 
reported  here,  if  defect  concentration  is  properly  taken  into  account. 

The  (visible)  optical  absorption  spectrum  of  these  samples  is  consistent 
with  a  charge  state  of  3+  for  the  Fe  .  Furthermore,  the  pressures  and 
temperature  above  apparently  result  in  a  Fermi  level  significantly  above 
the  Fe  *  trapping  level  (but  below  the  Fe2+  level) ,  since  the  <11%  concen¬ 
tration  depends  only  very  weakly  on  P}1  Q  (see  Ref.  1).  Thus,  since  the  Fe 
concentration  is  much  higher  than  that  of  the  accidental  impurities,  and 
presumably  there  are  no  other  donors  present  in  significant  concentrations, 
the  <H  >  and/or  <D  >  concentration  should  be  approximately  equal  to  the  Fe 
concentration.  Repeated  cycling  between  approximately  95%  <D%,  5%  <H%  and 
5“  <D  >*  95%  <H  *  (monitored  by  IR  absorption),  during  which  time  the  total 
<U  >  ♦  <l)%  concentration  remained  constant  with  experimental  error  (±  2%), 
indicated  a  total  <H  >♦<!)>  (and  hence,  Fe)  concentration  of  5.6  x  lO^/cm3. 

Since  the  IR  absorption  peak  shape  varies  considerably,  depending  on 

purity,  temperature,  etc.,  we  have  analyzed  our  data  in  terms  of  the  integrated 

absorption  (although,  except  for  calibration,  actual  measurements  were  made 

only  at  fixed  wavelengths,  as  discussed  in  Sec.  IV),  which  should  be  essentially 

independent  of  line  shape  and  hence  can  be  used  as  a  measure  of  <H%  or  <D% 

concentration  in  any  rutile  crystal,  independent  of  purity.  Let  D  (v)  be 

H 

the  optical  density  due  to  <H  >  in  a  particular  specimen  at  wavenumber  v 
and  let  be  the  corresponding  integrated  absorption,  i.e., 

A„  =  S  DH(v)dv  (1) 
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D(v)  is  related  to  the  absorption  coefficient  a(v)  by  the  relation 
D(v)  ■  a(v)t/2.3 

where  t  is  the  sample  thickness  in  cm.  Also,  we  define  a  quantity  A  j  (the 

"absorption  strength"  per  ion),  which  we  later  relate  to  the  oscillator 
strength,  as 


a 

H 


1 _ 

2 . 3n^ 


I  a(v)d\> 


(2) 


where  nH  is  the  <H*>  concentration  in  cm'3  (note  that  this  is  not  the  same 
notation  used  in  Ref.  1).  Further,  we  let  be  the  mass  of  a  proton,  MH 
the  total  mass  of  <H+>  in  the  crystal  and  H  the  total  mass  of  the  crystal. 


Then 


"h  ■  Vh  v 


(J) 


where  V  is  the  sample  volume  and  p  ■  4.26  gms/cmJ  is  the  density  of  rutile. 
Combining  Eq.  (1)  -  (3),  and  equivalent  equations  and  definitions  for  <D*>, 
and  letting  A  denote  changes  in  these  quantities  between  successive  runs, 
we  obtain 


M 


“  *  “h  4  “d  ’  o  (4"h  "h  *  4"d  "d> 


M  fMH  *H  ^2  % 

P  V  V 


(4) 


SJ<*h***V  * 


assuming  that  all  mass  change  was  due  only  to  <H*>  and  <D+>,  .is  was  the 
case  for  our  measurements.  K  contains  only  known  constants  and  geometrical 
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factors;  the  constant  8  (which  would  equal  2  if  the  absorption  strengths  of 
<H  >  and  <D+>  were  e<*ual)  can  be  determined  independently  of  the  mass 
measurements  if  chemical  equilibrium  is  maintained,  since  the  total  con¬ 
centration  of  <H+>  ♦  <D*>  will  remain  constant.  (In  principle,  some  departure 
from  this  condition  might  be  expected,  since  the  diffusion  coefficients  for 
<H*>  and  <D*>  should  differ  slightly,  but  in  practice,  this  difference  was 
not  observed,  presumably  because  changes  in  chemical  potential  which  would 
result  from  changing  concentration  constitute  a  much  stronger  driving  force 
for  diffusion  than  a  concentration  gradient  alone,  as  discussed  later  in 
this  section).  We  evaluated  8  from  absorption  measurements  taken  when  the 
crystal  contained  (essentially)  100%  <H*>,  and  100*  <D*>,  resulting  in  a 
value  for  8  of  3.78.  The  total  <H*>  ♦  <D*>  concentration  as  determined  in 
this  way  remained  constant  to  within  ±  2*  over  the  entire  range  of  <hW> 
concentrations.  Small  variations  were  to  be  expected  due  to  the  uncertainties 
in  partial  pressures;  hence,  we  feel  the  value  of  8  is  accurate  to  better 
than  ±  2*.  Assuming  that  the  uncertainty  in  the  determination  of  8  was 
smaller  than  the  uncertainty  in  the  determination  of  aH,  aH  was  then  evaluated 
from  a  least -squares  fit  of  a  series  of  6  successive  runs,  yielding  a  value 
of  (3.14  ±  .11)  x  10  ,  where  the  error  bars  represent  'J0%  confidence  limits. 

To  obtain  an  approximate  value  of  the  oscillator  strength,  f,  we  applied 
the  Smakula  equation3,  in  the  form 

f  -  2,3AH(D)MH(D)c2  On 

rH(D)  ,  — -  - — - 

716  1  nH(D)  C"2  *  2) 2 

wnere  e  is  the  electron  charge  in  esu,  n  is  the  refractive  index  at  the 
stretch  frequency,  c  is  the  velocity  of  light,  pH([))  is  the  reduced  mass  of 
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HfDI 

the  oscillator  and - - — *— < ^  is  the  area  under  the  absorption  coefficient 

_2 

peak  in  cm  .  Using  the  reduced  mass  of  free  OH  and  Of)  molecules,  and  a 
value  of  the  refractive  index  n  =  2.5  (which  is  the  average  for  the  two 
polarizations),  results  in  identical  values  for  f)f  and  fp  of  4.66  x  10’2. 

This  may  be  compared  with  values  obtained  for  K  and  Rb  halides  of  ff  ~ 

5  x  10  and  halides  of  ~  5  x  10  Thus,  the  oscillator  strength  in 
rutile  is  considerably  larger  than  that  in  the  alkali  halides  and  that 
calculated  foTtfus^ree  OH  ion.  liowever,  the  stretch  frequency  is 
significantly  shifted  te  lower  energies  in  rutile,  so  one  might  expect  a 
considerable  host  dependent  shift  in  f  also.  The  exact  agreement  of  f 

H 

and  fp  is  probably  accidental.  It  should  be  noted  that  the  values  of  af| 

and  6  arc  undoubtedly  more  accurate  than  those  for  fH  and  fp,  due  to  ambiguity 

of  the  meaning  of  uH^  and  uncertainty  in  the  appropriate  value  of  n,  as 

well  as  approximations  inherent  in  the  theory.  Thus,  the  quoted  values  of 

aH  and  B  should  be  used  in  experimental  determination  of  nfJ  and  »ip. 

As  indicated  in  a  previous  publication2,  the  <ll+>  line  shape  in 

rutile  is  rather  complex  and  depends  on  other  impurities.  The  measured 

peak  width  at  half. maximum  for  the  Fe-doped  crystals  used  here  was  -  27  cm-1 

for  the  <H  >  absorption  at  35#C,  which  was  somewhat  larger  than  for  undoped 
2 

crystals.  Comparison  of  the  measured  integrated  absorption  for  these 
crystals  with  that  for  a  Lorcntzian  line  of  the  same  width  and  height  results 
in  a  measured  area  10%  larger  than  for  the  corresponding  Lorcntzian  peak. 

In  addition  to  <11  >-associated  absorption  peaks  previously  reported2  at 
'  3650  and  4350  cm  1 ,  we  have  also  observed  an  absorption  peak  at  2119  cm’1, 
the  amplitude  of  which  is  proportional  to  <!!*>  concentration  and  is  smaller 
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than  the  main  <H*>  peak  at  3279  cm'1  by  -  200X.  <D%-associatcd  absorption 
spectra  centered  at  ~  2800  cm  1  (very  broad)  and  3200  cm'1  were  also  observed. 
(The  latter  is  difficult  to  observe,  since  it  is  "buried"  under  the  <H*> 
absorption).  These  spectra  apparently  correspond  to  the  4350  and  3650  cm*1 
<11  >  spectru',  and  have  similar  absorption  strengths.  The  4350  cm'1  absorption 
was  previously  tentatively  identified  as  the  1st  harmonic  of  the  3279  cm'1 
<H  >  peak.  It  now  appears  that  this  is  incorrect,  since  the  ratio  of 
4350:3279  is  almost  identical  to  3200:2438,  which  would  be  very  improbable 
for  such  strongly  shifted  harmonic  absorptions.  Thus,  a  combination  mode 
absorption  for  all  four  of  the  higher  energy  satellite  absorptions  appears 
to  be  the  most  reasonable  explanation.  This  would  yield  an  approximate 
excitation  energy  for  one  of  the  transverse  <H*>  modes  corresponding  to 
4350-3279  ■  1071  cm  1.  This  is  slightly  more  than  1/2  the  energy  of  the 
2119  cm  absorption,  which  thus  could  be  the  1st  harmonic  of  this  transverse 
mode.  Unfortunately,  the  predicted  fundamental  at  -  1071  cm"1  falls  in  the 
region  of  strong  lattice  absorption  and  is  probably  not  observable.  Polarization 
and  temperature  dependence  studies  of  these  spectra  are  planned. 

It  is  worth  noting  that  the  <H*>-<D*>  substitution  technique  described 
above  provides  an  excellent  (perhaps  the  only)  means  of  obtaining  accurate 
measurements  of  <H%  diffusion  in  such  a  crystal.  Previous  attempts  have 
been  made  to  measure  <H  >  diffusion5,  but  agreement  among  these  measurements 
has  been  quite  poor.  It  now  seems  clear  that  the  difficulties  with  these 
experiments  stem  from  two  related  sources:  first,  the  driving  force  for 
diffusion  is  actually  the  gradient  in  the  chemical  potential;  only  in 
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special  cases  can  this  be  adequately  approximated  by  the  concentration 
gradient,  as  is  usually  assumed.  In  materials  like  rutile,  in  which  the 
chemical  potential  depends  strongly  on  ue,  which  in  turn  depends  strongly 
on  the  concentration  of  <H+>,  this  may  be  a  very  bad  approximation.  Secondly, 
as  we  have  shown1,  heat  treatment  in  an  atmosphere  of  H2,  which  has  been  used 
in  most  experiments,  results  in  introduction  of  Ti  interstitials,  which 
can  seriously  affect  the  results  of  the  measurements  in  a  variety  of  ways. 
While  the  binding  energies  required  for  detailed  calculations  are  not  yet 
available,  it  is  clear  from  the  analysis  of  Ref.  1,  that  heating  in  H,  with 

• 9 

no  H20  present,  should  result  in  a  very  high  concentration  of  <Ti4%  (see 
Oq.  (22),  Ref.  1).  We  have  verified  this  experimentally  by  heating  an 
undoped  crystal  to  800°C  in  10  Torr  of  very  dry  H r  This  resulted  in  a 
substantial  concentration  of  donor  interstitials,  as  indicated  by  the  optical 
absorption  and  electronic  conductivity  (~  10"4/«  cm),  but  essentially  no  H 
in  the  crystal  (<  1017/cm3,  as  measured  by  IR  absorption.)  This  is  in  marked 
contrast  to  the  results  obtained  from  heating  similar  crystals  in  an  atmosphere 
of  H20  and  02,  which  produces  large  concentrations  of  H  interstitials  and, 
under  proper  conditions,  virtually  eliminate  Ti  interstitials,  as  discussed 
above.  The  technique  of  <H+>  -  <D+>  substitution  described  above  eliminates 
both  of  these  problems,  since  chemical  equilibrium  is  maintained  at  all 
times.  Actually,  this  is  not  strictly  correct,  since  the  diffusion  coefficients 
for  <H*>  and  <D+>  should  not  be  exactly  the  same;  our  experimental  data, 
however,  indicates  negligible  change  inthe<H+>  r  <D*>  concentration.  This 
observation  clearly  demonstrates  the  importance  of  contributions  to  the 
driving  force  for  diffusion  other  than  concentration  gradients,  since 
expected  departures  from  chemical  equilibrium  due  to  differences  in  diffusion 
coefficients  are  not  observed.  Thus,  the  electric  fields  produced  by  small 
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changes  in  ye  enhance  <D+>  diffusion  and  retard  <H+>  diffusion,  so  that  only 
the  average  of  the  two  diffusion  coefficients  can  be  observed.  Diffusion 
measurements  using  this  technique  are  in  progress  and  will  be  reported  in 
a  future  publication. 

The  authors  wish  to  thank  Drs.  John  Shaner  and  Fritz  Lut y  for  numerous 
helpful  discussions. 
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THE  BAND  GAP  ABSORPTION  OF  RUTILE  (Ti02)* 


John  W.  Shaner 
Department  of  Physics 
University  of  Utah 
Salt  Lake  City,  Utah  84112 


The  existing  optical  and  phonon  data  on  rutile  have 
been  examined  to  determine  the  band  gap  of  rutile  (TiO  ). 
Arguments  are  presented  for  the  conclusion  that  the  band 
gap  is  indirect,  involving  electronic  transitions  between 
r  and  M  or  Z  and  A,  and  that  the 


gap  energy  is  3.02  eV. 
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Any  attempt  to  calculate  the  band  structure  of  rutile  (Tit^)  has 

in  the  past  been  met  with  a  large  number  of  conflicting  statements 

concerning  the  gap  energies.  Early  optical  absorption  measurements 

showed  the  onset  of  Urbach  tail  absorption  starting  near  3  eV*,  and 

this  has  been  taken  as  the  gap  energy  by  many  authors.  Subsequent 

reflectivity  measurements  have  shown  that  a  significant  increase  in  the 

absorption  coefficient  begins  to  appear  only  at  energies  greater  than 
2 

3.5  eV  .  Interband  Faraday  rotation  and  two  photon  absorption  measure¬ 
ments  have  led  to  suggestions  of  an  indirect  gap  of  3.0  eV  with  a  direct 

3  4 

gap  of  3.3  eV  and  a  direct  gap  at  T  of  3.75  eV  ,  respectively. 

Electromodulation  spectra  have  been  measured  for  rutile^’^  with 
the  conclusion  that  there  is  optical  absorption  structure  around  3.05  eV 
which  is  consistent  with  phonon  assisted  transitions.  From  the  energy 
difference  between  the  phonon  absorption  and  phonon  emission  peaks,  an 
estimate  can  be  made  of  the  phonon  energies  involved:  30  meV  for  light 
polarized  E  j^c-axis,  and  40  meV  for  light  polarized  E| |c-axis5.  Since 
the  energy  halfway  between  the  phonon  emission  and  phonon  absorption  peaks 
is  the  same  for  both  polarizations,  we  may  assume,  as  in  Ref.  5,  that  the 
initial  electronic  states  are  the  same,  and  final  electronic  states  are 
the  same,  for  both  polarizations. 

Recent  neutron  scattering  expei.  „*nts,  elucidating  the  phonon 

structure^,  allow  one  to  fit  the  electromodulation  data  and  to  establish 

the  possible  initial  and  final  electronic  states.  Since  both  the  initial 

states  and  the  final  states  are  the  same  for  both  polarizations,  then  by 

k-vector  conservation  the  phonons  must  come  from  the  same  point  of  the 

7 

Brillouin  zone.  From  the  measured  phonon  dispersion  curves  ,  only  one 
such  point,  the  M  point  (zone  boundary  in  the  [110]  direction  -- 
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has  a  significant  phonon  density  of  states.  At  the  M  point,  there  are  two 
phonon  modes^with  the  symmetry8  or  M(.  degenerate  at  39  meV  and  a  phonon 
mode  with  M9*  symmetry  at  32  meV.  These  energies  agree  with  the  electro- 
absorption  results  to  within  the  quoted  experimental  uncertainty.  We  may, 
therefore,  assume  that  the  indirect  transition  near  3.05  eV  involves  an 
emitted  M  point  phonon  and  an  energy  gap  between  electronic  states  of  3.02  eV. 
K-vector  conservation  restricts  the  electronic  transitions  involving  initial 
and  final  states  at  high  symmetry  points  to  1)  r  -*■  M,  2)  M  r,  3)  X  X 
4)  Z  -*■  A,  5)  A  -*■  Z,  and  6)  R  -*■  R. 

The  number  of  possible  assignments  may  be  significantly  reduced  by 
calculating  the  allowed  phonon  assisted  transitions  between  the  various 
electronic  states.  The  correct  assignment  must  involve  an  Mg  or  phonon 
for  the  E)  |c  polarization  and  an  m/  phonon  for  the  E  1  c  polarization,  and 
it  must  have  the  same  initial  electron  state  and  the  same  final  electron 
state  for  both  polarizations.  Following  the  selection  rule  calculating 
procedure  of  Lax  and  Hopfield9,  and  the  character  tables  of  Gay,  et  al.8, 
we  have  calculated  the  selection  rules  which  satisfy  these  criteria10.  The 
X  X  transitions  and  the  R  -*■  R  transitions  can  easily  be  eliminated.  The 
allowed  indirect  X.  -  X.  transitions  involving  the  appropriate  phonon 
symmetries  and  light  polarizations  are  Xj  -  Xj  and  X2  -  x  .  Since  both 
of  these  are  also  allowed  direct  transitions  for  X  polarized  light,  one 
would  expect  a  strong  zero  phonon  absorption  in  addition  to  the  phonon- 
assisted  transitions.  Since  this  behavior  is  not  observed,  the  pertinent 
electron  states  are  not  at  X.  The  allowed  R  -  R  indirect  transitions  would 
not  be  allowed  direct  transitions,  but  the  M9+  phonon  assisted  transition 
would  be  allowed  for  both  polarizations  for  given  initial  and  final  states. 
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Thus,  the  polarization  dependence  of  the  electromodulation  data  rules 

out  R  ^  R  transitions.  The  remaining  allowed  transitions  subject  to  our 

restrictions  are  presented  in  Table  I.  If  one  assumes  that  the  valence 

bands  are  made  up  of  02  2p  states  and  the  conduction  bands  are  made  up 
.  4  + 

of  Ti  3d  states,  the  allowed  initial  and  final  states  are  further 
restricted  to  those  bands  made  up  from  these  atomic  states.  In  particular, 
this  selection  eliminates  rs+,  ,  Tf ,  and  T/  as  possible  final 

states  and  T 2  and  as  possible  initial  states11. 


A  suggestion  has  been  made  that  the  highest  valence  band  and  the 
lowest  conduction  band  at  T  have  symmetry4.  If  it  is  assumed  in 

addition  that,  as  in  SrTi0312,  the  highest  valence  and  lowest  conduction 
states  at  the  high  symmetry  points  are  on  the  same  sheets  as  those  at  F, 


then  the  compatibility  relations  require  the  M  states  to  have  M  M  or 
+  3*4’ 

Mg  symmetry,  and  the  Z  states  to  have  symmetry.  From  Table  I,  it  is 
clear  that  in  this  case  the  3.02  eV  electronic  transition  would  be 


M3,4’  M3,4 


r  +  z 

1 3  *  3 


A2,4*  °r  A2,4  +  Z3' 


Apart  from  the  symmetries  involved  in  the  indirect  transition,  there 


is  Svill  a  question  of  excitonic  structure.  Since  the  selection  rules  for 
transitions  to  an  S-exciton  are  the  same  as  those  to  the  continuum,  one 
might  expect  to  see  transitions  to  a  bound  excitonic  state  as  much  as 
several  tenths  of  an  eV  below  the  continuum.  In  fact,  additional  structure 
in  the  electromodulation  spectra  is  observed  at  3.10  eV5  and  at  3.3  eV6. 

If  the  3.02  eV  electronic  transition  were  to  an  S-excitonic  state,  and  one 
of  the  higher  transitions  were  to  the  continuum,  the  same  polarization 
dependence  of  the  energies  should  be  observed  for  both  transitions.  The 
3.10  eV  transition  energy  is  independent  of  polarization,  and  the  3.3  eV 
transition  is  clearly  observed  only  for  the  E | |  c-axis  polarization.  Since 
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neither  of  these  polarization  dependences  is  expected  if  these  transitions 
are  to  continuum  states  above  a  bound  excitun  at  3.02  eV,  we  conclude  that 
there  is  no  evidence  for  bound  excitons. 

These  considerations  lead  to  the  conclusion  that  the  structure  observed 
near  3.05  eV  in  electromodulation  spectra  of  rutile  arises  from  an  indirect 
electronic  transition  between  T  and  M  or  Z  and  A,  with  phonon  emission. 

The  energy  difference  between  electronic  states  is  near  3.02  eV.  This 
conclusion  is  consistent  with  electromodulation5 interband  Faraday 

3  i 

rotation  ,  and  band  edge  optical  absorption  data.  Since  the  electro¬ 
modulation  data  is  obtained  optically  from  an  intrinsic  crystal,  the 
renormalization  of  conduction  electron  energy  due  to  polaron  effects 
should  not  contribute  to  the  3.02  eV  energy  gap.  This  gap  will  then  be 
appropriate  for  the  rigid  ion  band  calculations.  Such  detailed  symmetry 
information  for  interband  transitions  in  rutile  has  not  previously  been 
available.  We  suggest  that  it  should  be  useful  as  a  constraint  on  any 
band  structure  calculation  for  rutile. 
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Table  I  Allowed  phonon  assisted  electronic 

transitions  involving  an  M.  or  M 

5  6 

phonon  for  £  polarized  light  and 
Mg  phonon  for  £  polarized  light. 


Initial  (final  State 

Final  (initial)  State 

V'  r4* 

_ 

V  M2 

r  +  r  + 

2  *  3 

m3.  m4 

r  “  r  ~  ; 

V  M6  j 

r  r  " 

2  »  *3 

M7'  M8  I 

r5+  i 

<  ! 

V  Z3 

: 

A4 

1 

Z2'  Z4 

;  V 

- f 

A3 

_  i 

_  7 
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FIGURE  CAPTION 


1.  Symmetry  points  of  the  Brillouin  zone  for  the  rutile  structure 
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"Vacuum  reduction"  of  rutile  in  the  temperature 
range  normally  employed  is  shown  to  be  caused  by 
residual  vacuum  system  impurities. 
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The  process  of  heating  rutile  (TIO7)  in  vacuum  to  achieve  non¬ 
stoichiometry,  and  the  associated  n-type  conductivity  and  blue  colora¬ 
tion,  has  been  used  by  scores  of  authors  in  many  different  investiga¬ 
tions  ,  and  is  one  of  the  "stand;  rd"  experimental  techniques  in  rutile 
research.  It  is  the  purpose  of  this  letter  to  show  that  this  so-called 
"vacuum  reduction"  is  in  fact  a  poorly  defined  and  uncontrolled  process 
resulting  from  contaminants  in  the  vacuum  system  (i.e.,  pump  oil,  water 
vapor,  etc.).  As  a  "dirt"  effect,  which  could  involve  impurity  doping 
of  the  crystal,  "vacuum  reduction"  is  obviously  not  a  desirable  exper¬ 
imental  technique.  Ln  fact,  the  term  "vacuum  reduction"  in  this  context 
is  quite  misleading,  since  the  process  does  not  occur  in  a  clean  vacuum, 
and  may  not  even  involve  "reduction",  which  implies  oxygen-def icient 
non-stoichiometry.  Since  the  term  is  in  common  usage,  however,  we  will 
continue  to  use  it  (with  quotation  marks),  with  the  understanding  that 
it  refers  to  whatever  processes  leue  to  n-type  conductivity  and  blue 
coloration  when  rutile  is  heated  in  a  dirty  vacuum,  and  that  it  may 
involve  doping  with  donor  impurities  rather  than  true  reduction. 

Qualitative  observations  that  both  the  rate  and  degree  of  "vacuum 
reduction"  were  strongly  dependent  on  the  quality  of  the  vacuum  system 
has  led  many,  though  by  no  metns  all,  of  those  studying  this  material 
to  the  view  that  "vacuum  reduction",  at  least  below  1000 °C,  is  strictly 
a  dirt  effect^.  To  the  best  of  our  knowledge,  however,  this  assumption 
has  not  previously  been  quantitatively  tested.  It  seemed  likely  that 
"vacuum  reduction"  was  not  qualitatively  different  from  treatment  in  an 
H2  atmosphere,  which  is  known  to  involve  both  H-doping  and  Introduction 


of  T1  Interstitials3.  Recant  work  o»  Cr34  -  associated  defect.  In 
rutile  ,  however,  appears  to  require  a  qualitative  difference  between 
these  two  processea,  and  stimulated  this  more  careful  Investigation 
of  the  nature  of  "vacuum  reduction". 

The  experiment  described  below  was  sufficiently  clean  to  estab¬ 
lish  that  most  of  the  reports  of  "vacuum  reduction"  have  resulted  from 
vacuum  system  contaminants.  A  crystal  of  nominally  pure  rutile5  was 
first  oxidized  in  pure,  very  dry  oxygen  at.  800°C.  The  crystal  was  then 
placed  in  an  O-ring-free  glass  system,  consisting  of  a  charcoal  trap, 

2  H sec  Vacion  pump,  and  fused  quartz  sample  tube,  '.he  major  aim  of 
the  system  design  was  to  avoid  pump  oil  contamination.  The  system  was 
rough-pumped,  then  sealed,  after  which  t’-.e  charcoal  trap  was  cooled  to 
77°K  and  the  rest  of  the  system  was  baked  at  200-300°C  for  three  hours. 
This  modest  bake-out  was  the  weakest  feature  as  far  as  cleanliness  is 
concerned.  After  bake-out ,  the  Vacion  pump  and  quartz  tube  were  seeled 
off  from  the  charcoal  trap  und  removed  from  the  bake-out  furnace.  In 

less  than  1/2  hour  the  Vacion  pump  had  reduced  the  system  pressure  to 
about  10~7  Torr. 


The  peak  optical  absorption  coefficient  at  1.55u,  F  c-axis,  was 
monitored  as  a  measure  of  the  degree  of  "reduction".  This  absorption, 
resulting  in  the  characteristic  blue  colorut ion  previous iy  referred  to, 
is  due  to  conduction  electrons  and  is  directly  related  to  tin  number  of 
defects  or  donors  introduced  into  the  crystal3.  The  optical  measurements 
weie  made  through  the  quartz  tube  at  room  temperature  without  breaking 
vacuum,  Figure  l.  shows  the  peak  opt  leal  absorption  coefficient , 
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temperature  of  the  sample  heating  furnace,  and  vacuum  in  the  sample  tube 
(measured  with  the  Vacion  pump)  as  a  function  of  sample  treatment  time. 
The  effects  of  out-gassing  of  the  quartz  tube  and  sample  are  clearly 
evident  in  the  vacuum  readings  and  in  the  slight  initial  coloration  of 
the  sample.  The  last  heat  treatment  was  done  in  a  different  furnace 
which  heated  a  larger  portion  ol  the  sampte  tube,  resulting  in  strong 
outgassing  from  a  portion  of  the  tube  not  previously  heated.  The 
optical  absorption  coefficient  at  1.55m  rose  to  over  loOenT1 .  The 
increase  in  pressure  precludes  the  certain  identification  of  the  strong 
coloration  at  1150°C  as  being  true  vacuum  reduction.  However,  the  fact 
that  there  was  no  observable  change  in  sample  coloration  during  23  hours 
at  over  1000°C  clearly  demonstrates  the  absence  of  true  vacuum  reduction 
at  this  temperature  and  below.  For  comparison,  a  typical  treatment  which 
involved  heating  at  820-850°C  for  30  minutes  in  vacuum  of  i0“2  Torr  has 
been  reported6  to  yield  an  optical  absorption  coefficient  of  13cm"1. 

The  nature  of  neither  the  slight  initial  "reduction"  nor  the  strong 
"reduction"  at  1150° C  is  known  at  present,  although  it  is  clear  that 
at  least  the  initial  coloration  was  due  to  residual  contamination.  It 
was,  however,  established  that  neither  of  the  "reduction"  processes 
resulted  from  H-doping.  The  H  concentration,  determined  by  1R  absorp¬ 
tion3  at  3279cm'1,  remained  below  detectable  limits  throughout  the 
experiment . 

The  question  of  what  really  happens  during  so-called  "vacuum 
reduction"  is  still  open.  It  could  be  that  the  "dirt"  fills  a  purely 
catalytic  role  producing  either  Ti  interstitials  or  0  vacancies,  or  that 


impurity  doping  of  the  crystal  takes  place.  This  question  is  under 
further  investigation. 

The  nature  of  the  donor  defects  introduced  at  1150°C  was  examined 

briefly  by  ESR.  After  heating  to  1150  C,  the  sample  yielded  an  intense 

ESR  spectrum  at  4°K,  similar  to  that  reported  by  Chester7  for  "vacuum 

reduced"  material.  An  attempt  was  made  to  oxidize  the  crystal  step  by 

.step  and  check  the  ESR  at  each  stage.  Unexpectedly,  a  25  hour  treatment 

in  100  Torr  of  oxygen  at  the  low  temperature  oe  490°C  removed  practically 

all  (>  90%)  of  the  coloration.  The  ESP  spectrum,  however,  was  practically 

unchanged.  The  remaining  coloration  was  quite  resistant  to  oxidation.  The 

oxidation  behavior  was  comparable  to  that  of  samples  containing  Ti  inter- 

34- 

stitials ,  but  no  Ti  ESR  spectrum  was  seen.  All  of  the  known  dopants  which 
diffuse  appreciably  at  490°C  (Li,  H,  Na)  can  be  eliminated  as  the  source  of 
the  easily  oxidized  component  of  the  coloration.  It  would  be  reasonable  to 
assume  that  a  defect  species  produced  by  true  vacuum  reduction  has  been 
observed.  It  is  possible  that  this  defect  is  involved  in  Lange's^  experiment 
on  microwave  acoustic  relaxation  in  rutile.  Lange  obtains  "vacuum  reduction" 
in  a  good  (10  7  Torr),  but  perhaps  not  clean,  vacuum  at  both  850°C  and  1050°C, 
but  only  observes  appreciable  acoustic  attenuation  after  the  higher  temperature 
treatment. 

On  the  basis  of  the  results  reported  here  and  other  work  on  the  thermo¬ 
dynamics  of  H  in  rutile  to  be  published  shortly,  it  is  strongly  recommended 
that  vacuum  reduction  be  avoided  as  a  technique  for  inducing  n— type  conduc¬ 
tivity.  Departure  from  stoichiometry  can  readily  be  achieved  by  diffusion 
of  evaporated  Ti  in  a  very  clean  vacuum  at  temperatures  of  -  800°C.  Alter¬ 
natively,  n-type  conductivity  can  be  induced  by  H-doping  at  lower  temperatures. 
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The  thermodynamics  of  this  process,  however,  is  quite  complicated,  and  can 
under  certain  circumstances,  result  in  production  of  Ti  interstitials  as 
well  as  H  interstitials.  In  most  cases,  if  the  doping  is  carried  out  in  a 
clean  ambient,  containing  both  H2  and  H.,0,  only  H  interstitials  result. 
This  problem  will  be  dealt  with  in  detail  in  a  forthcoming  publication9. 

The  authors  acknowledge  with  pleasure  the  useful  discussions 
with  their  colleagues,  in  particular  with  Dr.  John  Shaner. 
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Fig.  1.  Plot  of  peak  optical  absorption  coefficient,  sample 
heating  furnace  temperature  and  sample  tube  vacuum  vs.  sample  treatment 
time.  The  optical  absorption  is  for  light  of  1.55p  wavelength  with 
electric  vector  perpendicular  to  the  c-axis.  Each  dot  corresponds  to  a 
measurement.  The  lines  serve  only  to  connect  the  dots.  Vertical  segments 
indicate  readings  made  before  and  after  optical  measurements  at  room 


temperature. 
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